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ABSTRACT 


The common weathering sequence of feldspars at surface 
pre s s u re - temper at ure condition involves these steps, 
feldspar to mica, mica to clay and clay to gibbsite. On 
the otherhand, some recent reports suggest that a direct 
conversion of feldspar to gibbsite is possible bypassing 
the intermediate phases. This thesis reviews the physico- 
chemical conditions for direct g i bbs 1 1 1 z at 1 on of feldspars. 

+3 

Factors like structural changes, A1 Co- ordination, Eh-pH 
controls and thermodynamic viability have been discussed at 
length It was shown that direct conversion of feldspar to 
gibbsite requires special conditions like pH > 10, quick 
remo V a 1 of silica etc. 

Chemical mass balance of four soluble cations viz; Ha, 
K, Ca and Mg between parent rock, weathered product and 
ground water was used to calculate the time of weathering 
in several known bauxite deposits. It was found that the 
rate of weathering varies widely, from 2.3 metres per 
million years to as high as 1300 metres per million years, 
for same element in different rocks and also for different 
elements in the same rock,) 

A diffusion controlled computer model for weathering 
of feldspar to gibbsite was developed. Variables included 
ground water velocity, concentrations of H and dissolved 
silica species e.g., H 4 Si 04 and rock size. Effects of 
these parameters on the time required to attain steady 



state were evaluated For diffusion through a fi>ced 
distance of the order of 50 cms to 70 m , the time required 
to attain steady state was found to be 3,17 years to 77689 
years In comparison, the time of weathering of a rock 
having similar thicknesses calculated by mass balance was 
found to be of the order of few thousand years to few 
million years. This showed that diffusion plays a 
significant role only in the initial stage of weathering. 
After attainment of steady state it is likely to be 
controlled by the chemical reaction: 

-|.T 

KAIS 13 O 3 4 H + 7 H 2 O - AllOH)^ + K + 3H4S1O4 

amorphous 
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CHAPTER 1 


INTRODUCTION AND OBJECTIVES 

Weathering is a complex physico-chemical process of 
mechanical degradation and chemical alteration of primary 
rock forming minerals. Study of weathering processes is 
significant both from an academic as well as from an 
economic point of view. Many minerals of economic 
importance owe their origin to weathering processes. 
Notable among these is bauxite. 

Weathering reactions are, in general, extremely slow, 
incomplete and often irreversible (Krauskopf, 1979). One 
of the reasons for this slowness, is the pressure- 
temperature conditions at earth's surface (1 atm; pressure 
and 25°C temperature) under which most of the weathering 
reactions take place. 

It IS well known that bauxite deposits are dominated 
by three minerals, namely, gibbsite, diaspore and boehmite. 
Bauxite is a product of tropical weathering of feldspathic 
rocks. Existing views regarding genesis of bauxite can be 
broadly divided into two groups. 

According to one group, weathering of feldspars gives 
rise to gibbsite, but through an intermediate phase of 
clays (Valeton, 1972; Grubb, 1979; Kronberg et al., 1979). 
Another group advocates a direct derivation of gibbsite 
from feldspars (Mackenzie, 1958 quoted in Keller (1979), 
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Grubb, 1970; Keller, 1979), There are however, evidence of 
gibbsitization by both the processes in the s ame deposit. 
For example, Sen and Guha (1987) have reported direct 
gibbsitization of feldspathic rocks in older weathering 
profiles and desilication of kaolinite to gibbsite in 
recent prof lies. 

Opinion IS also divided on another aspect- According 
to some workers, weathering of feldspars involves removal 
of silica from the structure, leaving behind an Al-rich 
layer. Subsequent weathering rate depends upon the rate of 
diffusion of dissolved species through this layer. After a 
threshold value of silica concentration is reached in the 
external solution, back reaction is promoted, whereby 
kaolinite is formed (Wollast, 1967; Dobrovolsky, 1986). The 
other idea emphasizes total dissolution of feldspar in the 
associated solution and thereafter recombination of 
dissolved A1 and OH ions to give rise to gibbsite 
(Keller, 19 54; Krauskopf, 1979). A third mechanism favours 
selective dissolution of feldspar from structurally weak 
sites such as crystal edges, zone of dislocation etc. 
However, the Al-rich residue left locally on the feldspar 
surface is apparently too thin to be detectable by existing 
1 nstrument a 1 techniques (Berner, 1978; Berner and Holdren, 
1 978(a) and (b) ) . 

In view of the above lack of agreement on the 
mechanism of direct gibbsitization among previous workers, 
this thesis project was carried out with the following 
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objectives : 

1. To review the existing literature regarding direct 
g 1 bb s 1 t 1 z a t 1 0 n and g i bbs i t i z at i on through v n termed i ate clay 
phases 

2. To derive the precise physico-chemical constraints 
under which indirect and direct g i b bs i t i z at i on of feldspars 
can take place, from thermodynamic stability diagrams, 

3. To estimate the time of weathering of a given thickness 
of parent rock and its repl a cement by weathered product of 
equal thickness, by means of chemical mass balance of four 
soluble cations (viz» Na, K, Ca and Mg) and hence 
determining the rate of weathering. 

4. To utilize the boundary conditions derived above in a 
diffusion controlled weathering model of feldspar. This 
highlights the role played by rock size, percolation 
velocity of ground water etc on the time required to attain 
steady states for diffusion over fixed distance. 

5. To compare steady state time obtained by diffusion 
mechanism with calculated time of weathering and to draw 
relevant conclusions therefrom. 
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CHAPTER 2 

WEATHERING OF FELDSPARS 

2 1 WEATHERING PROCESS 

Most rock forming minerals are out of equilibrium at 
pressure and temperature conditions in which they are 
exposed at the earth's surface. Therefore, there is a 
natural tendency of tliese minerals to con\/ert to secondary 
minerals which are more stable in surface pressure and 
temperature conditions, i.e,, 1 atmosphere and 25° C 
respectively. This process of adjustment is collectively 
known as weathering. Leaving aside special contribution of 
biological processes, we can treat weathering as a chemical 
reaction between minerals and natural agents - principally 
water , 

The reaction of natural waters with rock forming 
minerals are further accentuated by the presence of 
dissolved CO 2 gas. A solution of CO 2 is mildly acidic. 

For example, it can be calculated that water saturated with 

—3 • 5 

atmospheric CO 2 at a partial pressure of 10 * atms* has a 

o 

pH value about 5.7 at 25 C (Garrels and Christ, 1965), 

This feebly acidic solution reacts with rock forming 

4 * 

minerals by a process of ion-exchange between H from 
solution and cations from the minerals. Thus Kao 1 i n 1 zat 1 on 
of potash Feldspar can be depicted as : 



aKAlSigOg + 2CO2 + IIH2O = Al2Si205(0H)4 + 2K*' 4- 41^43104 + 2HC 
2H'^ + 2HCO3 = 2CO2 + 2H2O 

2KAlSi30g + 2H'*' + WgO = = Al2Si205( OH )4 + 2K'^ + 4H4Si04 
On the basis of field observation and laboratory 


exper i men t s 

, common 

mafic and 

felsic minerals can be 

arranged in 

a series, 

accord i ng to 

their susceptibility to 

weathering 

( Kr auskopf , 

1 979 ) . 





Mafic 

Felsic 

Increasing 

i 


01 1 V 1 ne 

Ca^Na feldspar 

weather ab i 1 

ity 


Pyroxene 

Na-Ca feldspar 




Amph 1 bo 1 e 

K -f e 1 ds par 




B 1 0 1 1 1 e 

Muscovite 


Quartz 

It IS clear from the weather ab i 1 i ty series that in a felsic 
parent rock, feldspars are the more weatherable minerals. 
Among them, plagioclase feldspars are more weatherable than 
alkali feldspars. A c ommon weathering sequence is : 

Fe Idspars ^ Mica ^ Clay >> G i b bs 1 1 e 

Out of these four minerals, gibbsite 1s a hydroxide and the 
otherthree are silicates. 

In this thesis we are concerned with formation of 
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gibbsite by direct weathering of feldspars. It will be, 
therefore, necessary to discuss the weathering process in 
terms of crystal structure of the primary feldspars and the 
intermediate weathering sequence through mica and clay upto 
gibbsite. 

2.2 CRYSTAL STRUCTURE OF FELDSPAR 

Feldspars belong to the "Framework type" of silicates, 
in which S 1 O 4 tetrahedra are linked to one another by 
shared oxygens in all four corners, giving rise to a three 
dimensional network. For convenience of understanding the 
structure, the atomic arrangement may be conceived as 
linking of chains in two directions perpendicular to their 
length, although well defined chains of tetrahedra (as in 
pyroxenes and ampliboles) are not present (Deer, Howie and 
Zussman, 1979). The chains themselves are formed by the 
linking of inclined rings of four tetrahedra (Fig. 2. 1(a) 
and 2.1(b)). The repeat distance along the chain axis, 
i.e., a-crysta 1 1 ograph 1 c axis, is about four times the 
height of a tetrahedron. 

From a crysta 1 1 o~chemi ca 1 point of view, feldspars can 

be considered as modifications of quartz. When all four 

corners of SiO^ tetrahedron are shared, the formula of the 
4 

unit cell is Si (4, 1/2 0), i.e, S 1 O 2 In feldspars, 

+4 +3 

one-fourth of the Si positions are substituted by A1 

So, the formula becomes ^^ 3/4 ^2 AlSi^Og. 



A 


a 



° “ Oxygen ion 

• “ Silicon or Aluminium 
ion 


Fig. 2.1(a) Ideal illustration of feldspar chains 
(after Deer, Howie and Zussman, 1979) 



o -Oxygens 

MM-Mirror planes 


Fig 2.1(b) Part of sanidine (KAIS 13 O 0 , Monoclinic) 
structure viewed along normal to ( 001 ) 
(After Deer, Howie and Zussman, 1969) 
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+3 +4 

However, substitution of A1 for Si causes a 

deficiency of positive charge which is balanced by 
incorporating cations in the structure. Thus, in potash 
feldspars, the structural formula is KAIS 12 O 0 . Similarly, 
albite IS NaAlSi 303 Anorthite is CaAl 2 Si 203 where Ca +2 

balances the deficiency caused by substitution of two A1 

+4 

10 ns for two Si 10 ns Natural feldspars show solid 
solution series between these three end member composi- 
tions 

2.3 CRYSTAL STRUCTURE OF GIBBSITE 

The fundamental unit of the structure of gibbsite is a 

+ 3 

layer of A1 10 ns sandwiched between two sheets of close 


packed 

hydroxy 1 10 ns (Deer, 

Howie and Zu s sman , 

1 97 9 ) . 

Each 

A1 

lon IS octahedrally 

coord 1 nated , 1 . e 

, bonded 

with 

SIX neighbouring hydroxyl 

( OH ) 10 ns, three 

of which 

are 


situated above it and three below (part of this arrangement 

IS shown in Fig. 2. 2(a)). The gibbsite layer may be regarded 

as built up of octahedron linked laterally by sharing 

faces. In gibbsite, only two out of three octahedrally 

coordinated sites are occupied by A1 .Therefore, the 

+3 

structural formula is : A1 (6, 1/2 OH) i.e., AUOH)^ 0 ^ 

Al2(0H) ^ . A schematic diagram of gibbsite layer is shown 
in Fig. 2.2(b). 



n 



Fig. 2.2(a) Ideal structure 
X-axis 


@ Aluminium ions at 
©•Aluminium ions at 

QHydroxyl ions 

of gibbsite projected along 


Al2(OH)6 




Gi bbsite 


o Octahedral 
layer 


Fig. 2, 2(b) Schematic representat i on, of the structural 
unit of Gibbsite 
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2.4 CRYSTAL STRUCTURE OF KAOLINITE 

Clay minerals and micas are sheet silicates, where the 

basic structural unit is a silicate sheet formed by linking 

of three corners of the Si-0 tetrahedron leaving the fourth 

corner free. The structural formula of the building block 

IS therefore Si, (3, 1/2, 0), (1,0), i.e., (SiO or 

Si^O^* This unit has a net negative charge of >2, which is 

balanced in sheet silicates by a unique linking with 

+3 

modified gibbsite layer having Al in octahedral (i.e,, 

six-fold) coordination. For example, in Kaolinite 

+2 

tetrahedral sheet is linked to one (Al2(0H)^) octahedral 

layer in which two out of six OH corners of the gibbsite 

structure are replaced by oxygen atoms of the tetrahedral 
sheet (Fig. 2 *3 (a)). Thus, the structural formula of 

Kaolinite becomes Al2Sij^5(0H)^or AI2O3 .2Si0 2'2H^. Fig. 
2.3(b) shows the schematic representation of Kaolinite 

structure. 

2.5 CRYSTAL STRUCTURE OF MICA 

As stated above, micas are also s he e t - s 1 1 1 c ates with 

(Si^^ ) tetrahedral layers. However, in contrast to 

KaoliHhte, the building block is a three - 1 ayered T-O-T 

structure with one octahedral Al-0 layer linked on both 

sides with tetrahedral Si-O layers. As a result the basic 

“2 4*4 

formula of the Mica Group is - (S12O5) - (Al^ (OH) 2) 


Silicon 



© Aluminium 
O Oxygen 
© Hydroxyl 


Fig 2.3(a) Structure of kaolinite when Viewed along 

X-axis (after Deer, Howie and Zu s sman , 1 979 ) 


-2 

Si205 


Al2(OH)4 


Vw' 

0 


7-13 A 


Fig. 2.3(b) 


Schematic representation of the kaoHnite 
structure 
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-2 

(5i20 5^ ’ which actually represents the 

mineral pyrophyllite Micas have essentially the 

+ 4 +3 

pyrophyllite structure with substitution of Si by one A1 
ion. To restore charge neutrality, K**” ion is incorporated 

in the structure which occupies the interlayer position 

(Fig .2.4 ). The overall structural formula of muscovite thus 


can be wr i tten as 






IV 

SI3) 

^10 

(0H)2 

IV - Al'*’® 

in 

four 

fold 

co-ordination 

VI - Al'*'® 

in 

six 

fold 

co-ordination 


2.6 STRUCTURAL CHANGES DURING WEATHERING 

From the above descriptions of minerals, it is clear 


that the 

con ver s 1 on 

of 

feldspar to 

g 1 bbs 1 te 

through 

1 ntermed i ate 

steps of 

mica 

and Kaolinite is in 

fact 

a 

sequential 

format 1 on 

of 

successively 

simpler 

Si -0 

-Si 

linking in 

SI 1 1 c a t e s 

and 

ultimately 

/ 

formation 

of 

an 

octahedr a 1 1 y 

coordi nated 

A1 - 

OH structure. 

For e X amp 1 e i 



Four corner shared feldspar- — > three corner shared mica 
- — ^ three corner shared kaolinite — ^ No silica tetrahedra 
in g 1 bbs 1 te . 

+3 

There is a corresponding change in A1 co-ordination 

*^3 

also. In feldspars, all the A1 ions were in tetrahedral 
coordination. In mica (muscovite), two out of three A1 
ions in the structural unit are octahedrally coordinated. 
When the sheet structure of mica breaks down to give rise 


51205 ^ 

A1^0H)2 

5I26I 

51265 

aH0H)2 

S1265 



T 

7 


0 i 


/ 

T 



- -<K)- - 


\ 

T 

/ 

1 

0 


/ 

T 

k 


K - Potassium ion 


10 A° 


Fig. 2.4 


Schematic r epr - ^ en t a t 1 on of muscovite 
structure 
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to kaolioite, there is a further structural simplification. 

Two tetrahedral and one octahedrally coordinated layers are 

now reduced only to one tetrahedral and one octahedral 

+3 

layer. All the A1 ions present are now octahedrally 

coordinated. When kaolinite is weathered to yield 

gibbsite, the two-tier arrangement of octahedral and 

tetrahedral layers breaks to give way to a single 

octahedrally coordinated Al 2 ( 0 H )5 layer. 

The initiation of feldspar weathering takes place by 

the first surficial exchange of loosely bound K in the 

structure with H of the water. This renders the lattice 

weak and subsequent arrangement and adjustments take place 

in the structure, De Vore (1959) as quoted in Loughnan 

(1969), has discussed the mechanism of formation of 

1 ayered - 1 at 1 1 c e minerals, such as the micas, from feldspar. 

He pointed out that decomposition of the feldspar 

structure releases chains which have a certain degree of 

stability and retain the original Si-Al ordering of the 

tetrahedra. If these released chains are from the (100) 

and (010) surfaces of the feldspar, they can polymerize 

directly into tetrahedral sheets of composition (A), 

+3 

$ 13)0 ]_o Combination with octahedral cations (e.g., A1 , 
4*2 +2 “t3 + 

Mg , Fe and Fe ) and K would lead to the formation 
of clay minerals. However, if the released chains 
breakdown into individual tetrahedra, aluminium contained 
therein would be expected to assume its preferred 
octahedral coordination. In such cases, 1 ayered - 1 at 1 1 ce 
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+3 

silicates requiring at least part of the Al in the 
tetrahedral sheets could not develop 

It 15 pertinent, at this point, to mention that as we 
proceed from primary feldspar to secondary gibbsite via 
intermediates like muscovite and kaolinite, there is a 

successive decrease in Si;Al ratio. In alkali feldspars 
the ratio is 3 . 1 , in muscovite and kaolinite it is 1 1 , 

while in gibbsite it is zero. Thus, with progression of 

weathering, there is a progressive loss of not only cations 
but also silica from the structure. However, if in a 

particular case, Si Al ratio remains same in both parent 
and product mineral (e.g., in Anorthite ( CaAl 2 ^ 1 2*^8 ) 
kaolinite { Al 2S 1 2O OH ) 4} or muscovite { KA 1 3S 1 3O j_q( OH ) 2 to 
kaolinite) there is no silica loss from the parent mineral 
structure. In certain icases, where there is an increase in 
Si Al ratio from the parent mineral to the product mineral, 
silica has to be added in the system to ensure such 
transformations . 

For examp 1 e ; 

Anorthite —v Ca-montmori 1 Ionite 

CaAl2Si ^3 .17 2.33?’ 3.67 ^ 10 

(Raymahashay, 1984 ). These relationships are guided by the 

stoichiometry of weathering reactions in which Al is 

assumed to be immobile and is, therefore, locked up in the 

parent mineral. 

The above mentioned normal sequence of weathering of 
feldspars to kaolinite and kaolinite to gibbsite has been 
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recorded in numerous kaolinite and bauxite deposits 
(Valeton, 1972, Grubb, 1979, Kronberg, et al., 1979; Sen 
and Guha, 1987). However, weathering of feldspars directly 
to gibbsite is also recorded and is discussed further in 
the following paragraphs, 

2.7 EVIDENCES IN FAVOUR OF DIRECT GIBBSITE FORMATION 

In recent years, quite a number of field evidences 
brought to light the fact that direct g i bbs 1 1 1 z at 1 on of 
feldspar does take place in nature. Kaolinite, if present, 
IS only the product of re-silication of gibbsite. One of 
the earliest examples is that of Arkansas Bauxite Deposits, 
U.S.A., where frequent preservation of parent rock textures 
within the bauxite (gibbsite) was noted and conclusions 
were drawn in favour of direct alteration of feldspar to 
gibbsite (Mackenzie, 1958 as quoted in Keller (1979). 
Similar conclusion was also drawn by Grubb (1970) for the 
bauxites of Mitchell plateau, victoria, Australia, where 
basaltic (parent rock) textural features were preserved. 
Keller (1979) studied the Arkansas bauxite samples under 
scanning electron microscope and provided visual proofs of 
direct gibbsitic alteration of feldspars Very recently, 
Sen and Guha (1987) reported feldspar grains altered along 
the margins, directly to gibbsite, from the East Coast 
Bauxite Deposits, Orissa, India, 

In all the above cases, evidences of direct alteration 
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of feldspars to gibbsite, bypassing the intermediate 
kaolinite stage, are overwhelming. From a structural point 
of view, this IS somewhat unusual and enigmatic In the 
scheme of successive structural simplification through 
removal of silica and cations, this can be described as a 
"step-jump", where under certain combination of physico- 
chemical conditions, complex and tightly packed feldspar 
structure directly yields simple gibbsite structure. 

For this reason, direct g i bb s i t i z a t i on of feldspar can 
be considered as a special case. This thesis work is aimed 
at the details of special conditions which are needed to 
weather feldspars directly to gibbsite. 

2.8 KAOLINIZATION VERSUS G I BBS I T I Z AT ION - WITH 

+3 

REFERENCE TO A1 CO-ORDINATION. 

+4 

In contrast to Si in alkali silicate solution, 

which on neutralization changes from six to four fold co- 

+3 

ordination (Iler, 1955), A1 in alkali aluminate solution 

changes rapidly from four to six fold co-ordination on 

falling to pH 0. However, in presence of silicic acid 
+3 

(H 45104 ), four fold coordination may be stable down 

to pH3 (Grubb, 1970), Hence, the conclusion follows that 
at a pH value below 3, when six-fold A1 co-ordination Is 
fd'voured, gibbsite formation is accentuated with rapid 
removal of silica, When aluminium becomes mobile, it 
leaves the parent mineral structure to enter the solution 
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and eventually precipitates as A1(0H)2- Here, role of OH 

ions should be conspicuous in raising the ion activity 

+3 **• 

product of A1 and OH ions above the solubility product 

of A 1 { 0 H )3 (gibbsite), so that gibbsite could precipitate. 

According to Garrels and Christ (1965), gibbsite formation, 

unlike kaolinite formation is largely the result of 

+ 

unbuffered rise in pH (due to continuous uptake of H from 
the solution and removal of silica) during weathering. On 
the otherhand, if gibbsite were to be formed from 
kaolinite, it can simply form at any ordinary Intermediate 
pH by desilication and consequent structural rearrangement 
of kaolinite. 

It will be pertinent, at this point, to have a closer 

look on dependence of aluminium mobility on pH. fig. 2,5 

depicts the stability fields of various species of 

aluminium. The isoelectric point for the aluminium and 

aluminate ions is around pH6 On the left hand side of it 
+ 3 

AI is the predominant ion in solution, while on the 

right hand side of it, AUOH)^ is the predominant ion in 
solution. Furthermore, at around pH6, AMOH)^ (gibbsite) 
has the minimum solubility value. So, the effect of very 
high or very low pH is to mobilize aluminium effectively 
fr6m the parent mineral as soluble ions. Once aluminium is 
in solution, the solution pH should be centred around pH6 
in order to have maximum gibbsite precipitated. This 
sudden change in pH, however, cannot be accounted for with 
cer ta 1 nty . 



loq (concentration) 



Fig. 2 5 


Stability fields of various species of 
aluminium 
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All these evidences point towards the fact that with 

respect to pH range also, direct g i bb s 1 1 1 z a 1 1 on of 

feldspars is a special case, which is only possible at a 
very low or very high pH. 

2.9 KAOLINIZATION VERSUS G I B B S 1 T I Z AT I ON - WITH REFERENCE 
TO STABILITY DIAGRAMS 

Fig. 2.6 (drawn from the data of Appendix 1) shows the 
stability relationships among various phases of K-feldspar, 
kaolinite, muscovite mica and gibbsite. The diagram is, 
actually, a superposition of stability relationships among 
members of two systems. One is the system K-feldspar~ 

muscovite m i c a - k ao 1 i n 1 1 e - g i bbs i te , all being in crystalline 
state. The other is the system K-feldspar (crystalline) - 
muscovite mica (crystalline) - kaolinite (poorly crysta- 
lline)~gibbsite (amorphous). 

The various phase boundaries are as follows , 

1. Equilibrium boundary between crystalline kaolinite 
and crystalline gibbsite 
2 Equilibrium boundary between K-feldspar and 
crysta lime kaolinite 

3. Equilibrium boundary between muscovite mica and 
crysta lline kaolinite 

1'- Equilibrium boundary between poorly crystalline 
kaolinite and amorphous gibbsite 
2'- Equilibrium boundary between K-feldspar and 



Fig. 2 6 Stability diagram for Gibbslte. kaolinite, 

Muscovite and K-fe1dspar at 25°C 
' K-feldspar (crystalline) 

m - m 1 c a ( crys t a Mine) 
g - crystalline gibbsite 
- amorphous gibbsite 

k - crystalline kaolinite 

kp - poorly crystalline kaolinite 
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poorly crystalline kaolinite 
3'- Equilibrium boundary between mica and poorly 
crystalline kaolinite 

1(a) - Metastable extension of equilibrium boundary 1 

1(b) - Metastable extension of equilibrium boundary 1' 

2(a) - Metastable extension of equilibrium boundary 2 

2(b) - Metastable extension of equilibrium boundary 2' 

Lines 1(a) and 2(a) meet at point P and lines 1(b) and 
2(b) meet at point Q . 

4 - Metastable equilibrium boundary between K-feldspar 

and crystalline gibbsite 

5 - Metastable equilibrium boundary between K-feldspar 

and amorphous gibbsite. 

It can be seen in the diagram that in both the cases 
of crystalline and poorly crystalline or amorphous phases, 
K-feldspar and kaolinite have common boundary and so have 
kaolinite and gibbsite. This indicates that formation of 
kaolinite from K-feldspar and gibbsite from kaolinite is 
thermodynamically predictable This can occur under a wide 
range of pH variation along with silica removal. On the 
other hand, nowhere K-feldspar and gibbsite have a stable 
equilibrium boundary. This rules out the possibility of 
direct conversion of K-feldspar into gibbsite via stable 
equilibrium relationship. This, therefore, calls for 


metastable 

extens i ons 

of 

re levant 

phase 

bound ar 1 es , by 

which, we 

arrive at 

two 

points P 

and 

Q as described 

) 

earl i er 

Lines 4 and 

5 start from 

P and 

Q respect i ve 1 y. 
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So, P IS the lowest point on the line 4 which furnishes the 
composition of water, where K-feldspar and gibbsite (both 
in crystalline phases) are in nietastable equilibrium. On 
the otherhand, Q is the lowest point on the line 5 which 
furnishes the composition of water, where K-feldspar 
(crystalline) and amorphous gibbsite are in metastable 
equilibrium. 


or ^ ^ 


log 

■sr 

• 

CO 

SH'*' 

or, 

‘4 - 10=-“ 
10-= •* 

log 

ir+ = 6*2 


. p,6.2 

=10 , or, 


H 

H 

10"^ *2 




'K 


Taking the average concentration of in natural stream 

“4.23 

water as 2.3 ppm (= 10 m/1) (Mason and Moore, 1982), we 
get two different values of pH, viz; 12.63 and 10.43 
corresponding to P and Q respectively. Both indicate 
solutions of high alkalinity. This lends support to the 
contention that direct alteration of K-feldspar to gibbsite 
takes place at high pH conditions. 

In Na- and Ca- systems, however, direct conversion of 
albite and anorthite to gibbsite is thermodynamically 
stable. In Na-system (Fig. 2.7, drawn from the data of 
Appendix 1), the point R lies at the point of intersection 
of k ao 1 1 n 1 1 e - g 1 bb s 1 te and a 1 b 1 1 e -k ao 1 i n i t e equilibrium 
boundaries This is the lowest point on the equilibrium 
boundary between crystalline albite and crystalline 
gibbsite It furnishes composition of water, where 
crystalline albite and crystalline gibbsite are in stable 


g i bbsite 
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equilibrium Taking the average Na concen tr at i on in 
natural stream water as 6.3 ppm {= 10 m/1) (Mason and 
Moore, 1982), we get the pH of the water as follows 


From R, log = 10.20 




^Na'*’ ^ ,q10.20 


or, a^+ = X = j^q-13.76 


pH = 13.76 


In Ca-system (Fig. 2.8, drawn from the data of 

Appendix 1), the point S lies at the point of intersection 

of k ao 1 1 n 1 te - g 1 bbs 1 t e and anorth i te-k ao 1 i n i te (all in 

crystalline state) equilibrium boundaries. This is the 

lowest point on the equilibrium boundary between 

crystalline anorthite and crystalline gibbsite. It 

furnishes composition of water, where crystalline anorthite 

and crystalline gibbsite are in stable equilibrium. Taking 

the average Ca concentration in natural stream water as 
— 3 ♦ 42 

15 ppm (= 10 m/l) (Mason and Moore, 1982), we get the 

pH of the water as follows . 


From 

S, log 

= 17.27 or, 

^Ca+2 _ 

o ” 



<v> 

( a ^+)2 

or, 



= 10-20-69 

or, 

aH+ = 

j ^ O - 10.345 

= 10.345 


In both the above cases, 
highly alkaline range . 

These calculations indicate 


pH values deduced lie in 


that direct g i bbs i 1 1 zat i on 


log (a^^+2/(aHy) 


2A 



ig. 2.8 


Stability diagram for gibbsite, kaolinite 
and anorth i te at 25°C 
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of K, Na and Ca-feldspars requires highly alkaline solution 
pH value While Na and Ca feldspars can be in stable 

equilibrium with gibbsite, K-feldspar can convert to 
gibbsite only under metastable equilibrium conditions^ 

2.10 ROLE OF REDOX POTENTIAL (Eh) IN DIRECT G I 8B S I T I Z A T I ON 

Redox potential or "Eh" value of water has a well 

defined role in the formation of latesites, which are 

normally associated with bauxites. This is because, higher 

+2 +3 

Eh value helps in oxidizing Fe ions 'to Fe state, 

whereby it is precipitated in the form of Fe( 0 H )3 , and 

, mobility of Fe( 0 H )2 is much less compared to its ferrous 

counterpart. On the otherhand, any change in Eh of the 

+3 

soljUtion does not affect A1 ion and the process of 

I 

gibbsite formation. These aspects can be illustratec^ with 

the help of following diagrams (Fig. 2, 9(a) and (b)). F i g .2* 

9(a) shows the i r on - s t ab 1 1 i ty fields, where it can be 

+3 

observed that the field for ferric iron (Fe ) lies in the 
extremely high Eh range. On the contrary, the stability 
field of gibbsite covers a wide range of Eh from positive 
to negative (Fig .2.9 ( b ) ) . 

Furthermore, evidences in favour of gibbsite formation 
both under high and low Eh values can be found. For 
example, Keller et al (1954), as quoted in Grubb (1970), 
has found gibbsite formed under oxidizing conditions. On 
the otherhand, in the Gippland Bauxite Deposit, Australia, 


:h 0 

'o I ^s) 


-1 

0 7 14 

pH 

>"'9 2.9(a) Eh-pH diagram for the system Iron-water at 

25 C, after Garrels and Christ (1965) 

{ from Petersen, 1971) 



Fig. 2, 9(b) Eh.-pH diagram for the system Aluminium- 

water at 25^0, after Pourbaix, 1966 (from 
Petersen , 1971 ) 
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bauxite is in association with small amount of resicludl 
sulphide minerals, which indicates a reducing environ nient 
(Grubb, 1970). 

Besides the aforesaid chief factors, other factors 

like dense vegetation cover with intense rainfall 

(available in tropics), rugged topography and higher 

atmospheric temperature are also required to promote direct 

g 1 bb s 1 1 1 z a 1 1 on . Rugged topography helps in quick drainage 

of silica. Certain tropical forests can extract 

significant amount of silica from soil, thus help in silica 

removal (Lovering, 1959, quoted in Sen and Guha, 1987). 

According to Okainoto et al ( 1 957 ) as quoted in Sen and 

Guha, 1987, solubility of silica gel and the rate of its 

0 

dissolution increase three-fold from temperature 0 to 
0 

58 C Thus, higher ambient temperature also facilitates 
quick silica removal. Sen and Guha (1987) also emphasized 
the roles played by vegetation, rainfall and rugged 
topography in the direct g i bb s i t i z a t i on of feldspathic 
rocks in the East Coast Bauxite Deposit of India. 

To sum up, following points seem most important in 
considering direct g i bbs 1 1 i za t i on of feldspars as a special 
case as opposed to the existence of an intermediate 
kao 1 1 n 1 te stage : 

1. A sudden jump in the sequence of gradual simplification 
of compact feldspar structure 

2. Tendency of aluminium to change into six-fold co- 
ordination (from four-fold co-ordination in feldspars) at 
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very low (< 4) and at very high (> 10) pH solutions* 

3, Between pH values 4 and 10, aluminium practically 
remains immobile. Only at pH values outside this range, 
aluminium gains mobility in the form of water soluble ions 
and gibbsite can form from such a solution by precipitation 
at pH near the isoelectric point of aluminous species. 

4 Redox potential of water does not appear to play any 
conspicuous role in g i bbs 1 1 1 z at i on . 

5. Very high rainfall, good drainage and ample vegetation 
cover (as available in tropical regions) are needed in 
order to flush out silica and promote direct formation of 
gibbsite from parent feldspar. A fluctuating water-table 
can serve as a possible mode of silica removal. 

6. The following flow sheets furnish transformation 
relationships among different varieties of feldspare, mica, 
kaolinite and gibbsite, i.e., different pathways of 
gibbsite as well as kaolinite formation are shown. In all 
these transformations, at indicated pH, aluminimum is 
always considered to be immobile. 




Decreasing silica 


Resilication 


Ca-f eldspar 
( Anorthite) 



at pH > 10.345 
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CHAPTER 3 

TIME OF WEATHERING 

In the previous chapter, we have discussed the 
physico-chemical conditions for direct g i bb s 1 1 1 z a 1 1 on of 
feldspar. In the present chapter, we set out to 
investigate another quantitative aspect of feldspar 
weathering - that is - time. It is always interesting, 
both from an academic as well as from a practical point of 
view, to know the length of time involved to convert a 
given thickness of rock to its weathered product. Data 
from some known occurrences of b au x i te/ 1 at er 1 1 e are 
discussed below. 

3.1. Paragominas Bauxite Deposit, Brazil 

The feldspathic rocks have been weathered to bauxite 
in the Paragominas Bauxite Deposit, Brazil (Kronberg et al. 
1979). The time of weathering can be estimated from a 
cons j derat 1 on of chemical mass balance of four highly 
soluble cations, viz* Na, K, Ca and Mg. The relevant data 
for these elements are presented in Table 3 1 
Constraints : Average annual rainfall = 1000 m.m. = 1 m 
Problem : To find the time required to leach Na, K, Ca and 
Mg from a 2.3 m column of average crustal rock and convert 
into a 2.2 m thick laterite/bauxite profile. Here, the 
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TABLE 3.1 


E 1 emen t 

Average amount in 

Par agomi na s 

Bauxite Deposit 

( 

Average amount 

in earth's 

crust 

ng.g“^ 

Average amount 

1 n wor Id river 

water 

,(2) 

ng.g 

Na 


22700 

6 3 

K 

800 

18400 

2.3 

Ca 

800 

46600 

15 0 

Mg 

2500 

27640 

4.1 


Data from Kronberg et al (1979) 

(1) value after intense leaching 

(2) Data from Mason and Moore (1982) 
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thickness 2,2 m is the thickness of bauxite deposit in the 
Paragominas deposit. 


Procedure 


1 ) For Na • 

o 

Area of the base of the rock column = 1 m (say) 
Height of the rock column = 2.2 m 

3 

Therefore, volume of the rock column = 1x2.2 m 

= 2.2x10^ cm^ 


3 

Density = 2,7 gm/cm 

Therefore, weight of the rock column 

Na content in parent rock 
Na content in leached rock 
Therefore, amount of Na which 
contributes to water 


Therefore, total Na contributed by a 
= 5.94x10* x22550x10"° gin 

= 0. 134x10* gm 


= 2,7x2,2x10 gm 
= 5.94x10 gm 
= 22700 
= 150 

= (22700-150) 

= 22550 

= 22550x1 0"*^gm/gm of 
rock 

2.2 m column of rock 


-6 

In water, Na content = 6.3x10 gm/gm of water 
So , 


6.3x10 

Therefore 


” 6 > 

” gm of Na is in 1 gm of water 

, 0.134x10^ gm of Na is in 
12 

0.021x10 gm of water 
12 3 

= 0.021x10 cm of water 


In 1 year, total volume of water through rainfall 

2 3^ 

= 1mx1m =1m =10 cm^ 

Therefore, time required for weathering 

TO 

= 0.021x10 /10 years 

= 0.021 m . y . 


1 1 ) For Ic 

By the approach as before, weight of a 2.2 m 

6 

rock column = 5.94x10 gm 

"“d 

K content in the parent rock = 18400pg* g 

K content in the leached rock = 800^9* 

Therefore, amount of K which is contributed to water 

= ( 18400-800) p,g/gm of rock 

— 6 

= 17600x10 gm/gm of rock 

Total K contributed by a 2.2 m column of rock 

6 -6 
= 5.94x10 x 17600x10 gm 

= 0.105x10^ gm 

-6 

In water, K content = 2.3x10 gm/gm of water 

-6 

So, 2.3x10 gm of K is in 1 gm of water 
Therefore, 0.105x10 ^ gm of K is in 


12 


0.0456x10 

gm of water 

12 

o 

= 0,0456x10 

cm'^ of water 


In 1 year, total volume of water by rainfall 
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= 1 m X 1 = 1 = 10^ cm^ 

Therefore, time required for weathering 

= 0.0456x1 0^^ /10 ^ years 
= 0.0456 m , y . 

= 0 046 m y. 


Ill) For C a 

As before, weight of a 2.2 m column of parent rock 
= 5 . 94x1 0 gm 

-1 

Ca content in parent rock = 46600 M-cf • g 
Ca content in leached rock - 800 p-g.g ^ 

Therefore, Ca contributed to water 
= ( 4 6600 “ 800 ) i-tg /gm of rock 
= 45800x10 gm/gm of rock 

Therefore, total Ca contributed by 2.2 m rock column 
= 5.94x10^ x45800x1 o”^ gm 
= 0.272x10 gm 

In water, Ca content = 15l^g/gm of water 
So , 

“6 

15x10 gm of Ca is in 1 gm of water 

6 

Therefore, 0.272x10 gm of Ca is in 

12 

0.018x10 gm of water 
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12 o 

= 0 018x10 cni'^ of water 
In 1 year, total volume of water by rainfall 
== 1 =10^ cm^ 

Therefore, time required for weathering 

= 0,018x10^^ /10^ years 
= 0.018x10^ years 
= 0 » 018 m.y. 


iv) For Mg • 

As before, weight of a 2.2 m column of parent rock 
= 5.94x10 gm 

Mg content in parent rock = 2 7640 M-Q/gm of rock 
Mg content in the leached crust = 2 500 fig /gm of rock 
Amount of Mg contributed to water 
= ( 27640-2500) I^g/gm rock 

= 25140 M-^/gm of rock 

"• 6 

= 25140x10 gm/gm of rock 

Therefore, total Mg contributed by a 2.2 m high column of 
parent rock 

6 —6 
= 5.94x10 x25140x10 gm 

= 0.15x10° gm 

In water: Mg content = 4.lHg/gm of water 

= 4.1x1(r^ gm/gm of water 
So , 
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-6 

4.1x10 gm of Mg is in 1 gm of water 
6 

Therefore, 0 15x10 Qm of Mg is in 

( 

0.0366 X 1 O^^gm of water 
= 0.0366x10^^ cm^ of water 
In 1 year, total volume of water by rainfall 
= 1 =10^ cm^ 

Therefore, time required for weathering 
= 0.0366x10^^ /1C^ years 
= 0.0366x10^ years 
= 0 . 0366 m . y . 

= 0.037 m . y . 

The above calculations of time of weathering can be listed 
below , (Table 3.2) 

From the above calculations, it can be seen that the time 
required for leaching of a 2,2 m high column of crustal 

rock to give rise to a deposit of laterite/bauxite has a 

narrow range. In this case, this range lies between 0.018 
to 0.046 m.y. This, for all practical purposes, assigns 

the time required for weathering of a 2.2 m thick crustal 

rock at the Paragominas deposit, between 0.02 to 0.05 
million ye ars . 

With the availability of accurate field data, this 
simple technique can be used in other cases also, to have 
an ;dea of the time involved in the process of weathering. 
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TABLE 3.2 


For element Calculated time of weathering 


(m.y. ) 

Na 

0.021 

K 

0.046 

Ca 

0.018 

Kg 

0.037 
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3*2. East Coast Bauxite Deposit, India 

In the present section, we undertake another case 
study, where the relevant data are more accurate because 
composition of weathered rock and ground water are also 
available. The deposit in the "East Coast Bauxite 
Deposit", of Orissa, India. The particular areas of 
interest are Pottangi and Panchpatmali villages of Koraput 
district. Sen and Guha (1987) conducted a detailed study 
of the mature weathering profile developed over a 
khondalitic parent rock From surface, upto a depth of 
about 21 m, they delineated a zone - the so called 
"laterite-bauxite" zone, where according to them, evidences 
in favour of direct in situ g i bbs 1 1 1 zat i on of Khoudalite 
are plenty. From a depth of 21 m upto a depth of about 70 
m, the zone is called "saprolite" zone. This is also a 
zone of weathered/partially weathered material where 
evidence of directly altered gibbsiLe is less and amount of 
clay IS predominant. Below this zone, lies the unaltered 
parent rock (Khondalite), In the subsequent calculations, 
we take the thickness of weathered zone as 70 m, upto the 
contact between the saprolite zone and unaltered parent 
rock . 

Constraints : 

Average annual rainfall = 1.6 m = 1600 mm 

Depth of weathered zone below surface (i.e , depth taken 
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TABLE 3 3 


E 1 emen t 

Average 

Average 

Average amount 


amount in 

amount in 

in local ground 


parent rock 

bauxite bearing 

water from bau- 


( k hon d a 1 1 1 e ) 

weathered 

xite bearing 



prof 1 1 e 

hor 1 zon s 


( .g ) 


( tig.g ^ ) 

Na 

15525 74 

1990.75 

0.537 

K 

47646.50 

371.08 

0. 184 

Ca 

3466.00 

tr 

0.868 

Mg 

2133.81 

tr 

0.433 
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for calculations) = 70 m 

The data for the four elements are presented in Table 3,3. 

All the data are from Sen and Guha (1987). Except 
ground water analyses, all other values have been 
recalculated from percentage analyses. 

Procedure for determination of time of weathering • 

1 ) For Na 

Area of the base of the rock column = 1 m^ (say) 

Height of the rock column = 70 m 
Therefore, volume of the rock colum 

= 70m xlm^ =70m^ 

= 70 X 10^ cm^ 

Taking average density of the rock to be about 2.7 gm/cnr^ 

6 

weight of the rock column = 2.7x70x10 gm 
Now , 

“i 

Na content of the parent rock = 15525. 74 M'Q* g 
Na content of the leached rock = 1 990 . 7 5 • g ^ 

Therefore, amount of Na contributed to water 

= ( 1 5525 . 74-1990.7 5 ) ^ig. g“^ of rock 
= 1 3534. 99p-g /gm of rock 
= 1 3534 . 99x 1 0"**^ gm/gm of rock 
Therefore, Na contributed by a 70 m column of parent rock 
= 2.7x70x10^ X13534. 99x10“^ gm 
= 2.558x10^ gm 

In ground water, Na content = 0 537 fxg*g 

= 0.537x10" gm/gm of water 
So , 


I 
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0 537x 1 0 “^ gm of Na is in 1 gm of water 


Therefore, 2 558x10 gm of Na is in 


12 


4. 763x10 

gm of water 

12 

3 

4. 763x10 

cm of water 


In 1 year, total volume of water through rain fall 
= 1.6x1 rr^ = 1 6x10^ cm^ 

Therefore, time required i 

= 4. 763x1 0^^ /I .6x10^ years 
= 2.977x1 0^ years 
= 2.977 million years (m.y.) 

3.0 m.y. 

1 1 ) For K 

By the procedure shown as before, weight of a 70 m high 
column of khondalite = 2.7x70x10 gm 
Now , 

K content of the parent rock (khondalite) = 47646.50 
K content of the leached rock = 371. 08 
Therefore, amoun t of K contributed to water 
= (47646 50-371.08) ng.g“^ 

= 4727 5.42 Jig/gm of rock 
= 47275.42x10"^ /gm of rock 

Therefore, total contributed by a 70 m column of parent 

6 ~6 

rock = 2 7x70x10 x47275. 42x10 gm 

6 

= 8.935x10 gm 
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In ground water, K content = 0.184 

-6 

= 0 184x10 gm/gm of water 
So , 

-6 

0 184x10 gm of K is in 1 gm of water 

✓ 

Therefore, 8 955x10 ^ gm of K ' ' 

12 

48 56x10 gm of water 
= 48 . 56x1 0 ^ ^ cm ^ of water 
In 1 year, total volume of water through rainfall 
= 1 ,6x1 =1 .6x10° cm^ 

12 6 

Therefore, time required = 48.56x10 /I. 6x10 years 

= 30.35x10^ years 

f 

= 30.35 m. y . 

Ill) For Ca . 

As before, weight of the 70 m column of khondalite 

= 2.7x70x10^ gm 

Now , 

Ca content of the parent rock (khondalite) = 3466.00 
Ca content of the leached portion i.e., bauxite 

material = 0 (actually It was reported to be 
"trace" in amount, which can be neglected for the 
present calculation) 

Therefore, amount of Ca contributed to water = 3466 |ig .g 

—6 

= 3466x10 gm/gm of rock 
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Therefore, total amount of Ca contributed by 70 m column of 

6 

parent rock = 2 7x70x10 x3466. 00x10 gm 

= 0,655x10^ gm 

-1 

In water, Ca content= 0.868 1^0# g 

-6 

= 0.868x10 gm/gm of water 


So , 


-6 

0.868x10 gm of Ca is in 1 gm of water 


Therefore, 0.655x10 gm of Ca is in 


12 



0 . 7546x1 0 gm 

0 f 

water 

12 3 



0. 7546x1 0 cm 

of 

water 


In 1 yr., total volume of water by rainfall 

6 3 

= 1.6x1x10 cm 

. 6 3 

=1.6x10 cm 

12 6 

Therefore, time required = 0.7546x10 /1 . 6x10 years 

= 0.47x10^ years 
= 0,47 m . y . 

1 V ) For Mg 


As before, weight of a 70 m column of Khondalite 

‘ = 2.7x70x10^ gm 

Now , 

Mg content in parent rock (khondalite) = 2133.81 
Mg content in leached bauxite portion = 0 (reported to 
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be "trace" in amount) 

Therefore, amount of Mg contributed to water 

= 21 33 81 

-6 

= 2133 81x10 gm/gm of rock 

Therefore, total amount of Mg contributed by the 70 m 

6 

column of parent rock = 2.7x70x10 x2133, 81x10 gm 

= 0 4033x10^ gm 

In water, Mg content = 0.433 M-^-g ^ 

-6 

= 0,433x10 gm/gm of water 

-6 

So, 0.433x10 gm of Mg is in 1 gm of water 

Therefore, 0.4033x10^ gm of Mg is in 

12 

0 9314x10 gm of water 
12 S 

= 0.9314x10 cm of water 
Now, in 1 yr , total amount of water by rainfall 

6 3 

= 1.6x1x10 cm*^ 

12 6 

Therefore, time required = 0,9314x10 /I. 6x10 years 

6 

= 0.5821x10 years 
= 0. 5821 m.y . 

0.58 m.y. 

So, on the basis of geochemical mass balance of four commo n 
soluble cations, the time required to weather a 70 m high 
column of khondalitic rock and to replace it by a column of 


38 


TABLE 3.4 


For el emen t 

Calculated time 

(in million years) 

N a 

3.0 

K 

30.35 

Ca 

0.47 

Mg 

0.58 
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bauxite and clays of the same height can be listed below 
(Table 3.4) 

However, in all these estimates, we have considered that 
all the water furnished by rain has infiltrated below. 
This IS, however, not true as surface run off will 
constitute a large part of the total rain fall. Thus, 
effectively our age estimates are underestimates. The 
actual time of weathering will be more than what has been 
obtained These can be calculated if reliable data on 
percentage percolation are available 

DISCUSSION 

From Table 3.4, it can be seen that the age of 
weathering, for same thickness (70 m), calculated on the 
basis of four soluble catious vary among themselves. For 
Ca, it is least, being 0.47 m.y. while for K. it is 
largest, being equal to 30.35 m.y, so, the question that 
naturally arises, why is this discrepancy? 

As the measurements of concentrations were carried out 
by Sen and Guha (1987) with utmost caution, inaccuracy of 
data cannot be a possible answer. 

The explanation involves the aspect of relative 
mobility. Sen and Guha (1987) have measured the relative 
mobilities of all the catious by the ratio : 

/ cation in ground water 

Relative mobility = 7! , , . . TTtTT” 

•/ cation in average khondalite of tne 

area 
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TABLE 3.5 


E 1 emen t 

Relative mobility 

in ground water 

Time of weather i ng 

(m.y. ) 

Na 

3 495 

3.0 

K 

0.445 

30.35 

Ca 

20.00 

0.47 

Mg 

19.295 

0.58 
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The average relative mobilities for Na, K, Ca and Mg in 
ground water samples of Pottangi and Panchpatmali areas are 
presented below, along with the time of weathering 
furnished by each ion, in table 3.5 

From the above table, it can be seen that the relative 
mobilities are inversely proportional to the "time" 
furnished by corresponding ions. This probably holds the 
key to the solution of the problem of 1 1 me - d i s cr epancy . 

It IS well known that mobility of cations depend on 
two factors 

(a) the way a cation is bound up in a mineral mineral and 
the mineral stability, and 

(b) rate of precipitation after a cation is released from 
the parent mineral due to weathering. With these, another 
factor can be coupled in order to get the true magnitude of 
relative mobility. This is the retention of ions in the 
system. 

Three out of the above four ions, i.e., Ca, Mg and Na 
are released from their parent minerals more easily than K 
This is because, a f erro-magnes i an mineral or calcic 
plagioclase has higher weatherability than that of potash 
feldspar. Moreover, the present case is an open system, 
where ground waters are collected from flowing springs. 
Hence, it is plausible to argue that at a certain point of 
time of weathering, there are much more quantity of Ca, Mg 
and Na in ground water, than K. Again among Ca, Mg and Na, 
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Ca and Mg exceed the quantity of Na. This is evidenced by 
the values of their relative mobilities in water. However, 
as the system is an open one, more and more Ca, Mg and Na 
are irretrievably lost compared to K Thus we get lower 
" 1 1 me “ V a 1 ue s *' by Ca, Mg and Na but significantly higher 
value by K This contention leads to the fact that the 
concept of geochemical mass balance seems more applicable 
in case of the balance of K than the other three ions, in 

this case Hence, the time required to weather a 70 m 

column of khondalite seems to be around 30 m.y. 

With a different reasoning, Sen and Guha (1987) put 

the age of the 1 at e r i te - b au x i te duncrust as 50 + 5 m.y. 


The 1 r arguments 

are entirely 

qua 1 1 1 at 1 ve 

and 

speculative. 

Hence, the 

age 

determined on 

the basis 

0 f 

K 

s eems 

more 

acceptable 

as 

the real time of 

format i on 

of 

a 

7 0 m 

thick 


weathered prof lie. 

3. 3. Evidence from Microscopic view, East Coast Bauxite 

In this section, we will be applying the same 
principle as shown in Section 2 in order to calculate time 
of weathering in a much smaller scale. The principle of 
chemical mass balance is equally applicable in case of 
weathering of individual grains of feldspars. Figure 3.1 
shows a single crystal of feldspar which has been altered 
to gibbsite directly, along its margins (figure from Sen 
and Guha, 1987). Now, from the principle of chemical mass 
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TABLE 3,6 

For el erne n t 

Calculated time (years) 

Na 

4.63 

K 

46 82 

Ca 

0.725 

Mg 

0.9 
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balance it follows that other factors remaining constant, 
the length of time furnished by each element is 
proportional to the thickness of weathered layer. We can 
directly apply this principle in the estimation of time of 
weathering of a single feldspar crystal. Fig 3.1 shows the 
altered layer in thin section. For accuracy, this layer of 
gibbsite is measured in four places and the average 
thickness is taken for calculation 
Average thickness = 0 81 cm 

Magnification = 75 times 

Therefore, actual thickness= 0.81 cm/75 = 0.0108 cm 
The procedure to obtain the time of weathering is 
illustrated below for K only. For other ions procedure is 
exactly similar. 

In the previous section, we have obtained the time for 
weathering of a 70 m thick parent rock as 30.35 m.y., i.e., 
30.35x10 years on the basis of K-leaching, 

For weathering upto a thickness of 70 m, required time is 
30.35x10^ years 

For weathering upto a thickness 0.0108 cm, required time 


= 46,82 years. 

Table 3.6 shows the time furnished by each ion, obtained by 
the procedure shown above. 

Now, by the same arguments as have been outlined in Section 



UF - Unaltered feldspar 
Gibb^p- Cryptocrystalline gibbsite 


Fig. 3.1 Feldspar crystal replaced by cryptocrystalline 

gibbsite along margins {X75) (After Sen and 
Guha, 1987) 
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2, we can accept the time furnished by K 

So, the time required to create a 0 0108 cm thick band 
of gibbsite around the feldspar grain seems to be around 47 
years, in the East Coast Bauxite Deposit, Orissa, India, 
(Pottangi and Ponchpatmali plateau areas) This fact can 
be schematically represented by the diagram shown in Fig 
3.2. 

3.4 KARNATAKA BAUXITE DEPOSIT, INDIA 

Precambrian granitic gneisses (Tonalite gneisses) have 
extensively weathered to 1 ater i te-baux i te profiles, in the 
South Kanara district of Karnataka (Khanadali and Devaraju, 
1987) There are evidences of development of gibbsite, 
both by the direct alteration of feldspar of the parent 
rock and by desilication of kaolimte 

We have chosen the same four soluble cations, viz. Na, 
K, Ca and Mg for chemical mass balance. The quantities of 
these elements in parent rock and in weathered profile are 
recalculated from percentage analyses given by Khandali and 
Devaraju (1987). In the absence of ground water data, 
average world ' river water composition was used for the 
values of Na, K, Ca and Mg in water, from Mason and Moore 
( 1 982 ) . 

Constrai nt s : 

Average annual rainfall = 3500 mm (Ranganna et ai.^ 1984; 

Krishna Rao, 1971) 

= 3.5 m 


--t=0 



UF- Unallered 
feldspar 


Surrounding 
altered layer 
(gi bbsite) 


t=47years 


Fig. 3.2 



Idealised diagram of weathering of a feldspar 
grain showing time of weathering 
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TABLE 3.7 


Element Average amount 

in parent rock 

( Tonal 1 1 e 

gneiss) 

{ >9 ^ ) 

Average amount 

in g 1 bbs 1 te 

bean ng wea- 
thered prof 1 1 e 

) 

Average amount 

1 n wor 1 d r i.ver 

water 

, -1 , 

( *9 ) 

Na 

32633.06 

755.48 

6.3 

K 

12193.56 

1170.43 

2.3 

Ca 

19339.00 

205.13 

15.0 

Mg 

9429.80 

4890.95 

4.1 
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Depth of the weathered zone below surface (i e., depth 
taken for calculations) = 26 m 

The data for the four elements are presented below in 
tbaular form {Table 3.7), 

Procedure for determination of time of weathering . 

1 ) For Na 

p 

Area of the base of the rock column = 1 m (say) 

Height of the rock column = 26 m 

q 

Volume of the rock column = 26x1 m 

= 26x1 0^ cm ^ 

3 

Taking average density of the rock to be about 2.7 gm/cm 
Weight of the rock column = 26x10^ x2.7 gm 
Now , 

Na content of the parent rock = 32633.06 

Na content of the weathered rock = 755.48 lagig””^ 

Therefore, amount of Na contributed to 

water = (32633.06-755.48) M-g*g ^ 

-1 

= 31 87 7 58 M-g- g of rock 
= 31877.58x10 ^ gm/gm of rock 

Therefore, total Na contributed by a 26 m column of parent 
rock = 26x10^ x2 . 7x3 1 87 7 . 58x1 0"^ gm 
= 2.238x1 0 ^ gm 

In ave. world river water, Na content = 6.3 pg * g ^ 

-6 

= 6.3x10 gm/gm of water 

-• 6 

So, 6.3x10 gm of Na is in 1 gm of water 

Therefore, 2.238x10 gm of Na is in 
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0,355x10^^ gm of water 
= 0 355x10^^ cm ^ of water 
In 1 yr., total volume of water through rain fall 
= 3 5x1x1 O^cm^ 

Therefore, time required for weathering 

= 0 355x10^^ /3»5x10^ years 
= 0,1x10^ years 
= 0 1 m . y . 

1 1 ) For K . 

By the procedure shown in case of Na, weight of a 26 m high 
column of granite gneiss (parent rock) 

= 26x 1 0^ x2 . 7 gm 

Now , 

K content of the parent rock = 12193. 56 |ig*g ^ 

K content of the weathered rock - 1 170.43 

Therefore, amount of K contributed to water 

= (121 93.56-1 170.43) M-g. g 
= 11023 13 pg.g 

= 11023.13x10 ^ gm/gm of water 
Therefore, K contributed by a 26 m column of parent rock 

= 26x10°x2. 7x1 1023. 13x1 0~°gm 
= 0.774x10° gm 
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In ave world given water, K content = 23 

-6 

= 2.3x10 gm/gm of water 
So , 

— 6 

2,3x10 gm of K is in 1 gm of water 
Therefore , 

0.774x10 gm of K is in 

0.336x10^^ gm of water 
= 0. 336x10^ of water 

In 1 yr , total volume of water through rainfall 
= 3.5x1x10^ cm^ 

12 6 

Therefore, time required = 0.336x10 /3.5x10 years 

= 0 096x10^ years 
= 0.096 m . y . 


Ill) For Ca : 




As before, weight of 26 m column of parent rock 

= 26x 1 0 ^2 7 gm 


Now , 

Ca content of the parent rock = 1 9339.0 M-9 -g 

Ca content of the weathered rock = 205,1 3 M-g^g 
Therefore amount of Ca contributed to water 
= (19339.0-205.13) 

= 19133. 37|ig. g 
= 19133.87x10 gm/gm of rock 
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Therefore, total amount of Ca contributed by 26 m column of 
weathered gibbsitic product 


6 -6 
= 26x10 x2. 7x19133.87x10 gm 

= 1.343x10° gm 

In ave; world river water, Ca content = 15 M-g ♦ g 

-6 

= 15x10 gm/gm water 


So , 


-6 

15x10 gm of Ca is in 1 gm of water 


6 

Therefore 1,343x10 gm of Ca ' 


12 


0,0895x10 

gm of water 

12 

3 

= 0.0895x10 

cm of water 


In 1 yr, total volume of water by rainfall 


3,5x1x10 cm'^ 


J. z ^ 

Therefore, time required = 0,0895x10 /3.5x10^ years 

= 0.0256x10^ years 
= 0 . 0?56 m . y . 


1 V ) For Mg • 

As before, weight of a 26 m column of parent rock 
= 26x10^x2. 7 gm 

Now , 

Mg content in parent rock = 9525.8 lig • g 

Mg content in weathered rock = 4890.95 ^0*g ^ 
Therefore, amount of Mg contributed to water 


APR 1989 

CENL''<'\L L''=>RARY 

= (9529.8-4890. 

= 4638 .85 M-g *9 

— ft 

= 4638.85x10 gm/gm of rock 

Therefore, total amount of Mg contributed by the 26 m 

6 -6 

column of parent rock = 26x1 0 x2 . 7x46 38 . 85 x 1 0 gm 

^ 6 

= 0.326x10 gm 

-1 

In water, Mg content = 4.1 M-g-g 

“6 

= 4 1x10 gm/gm of water 

-6 

So, 4.1x10 gm of Mg is in 1 gm of water 


95) 

1 


-1 


I I 


K.\^ P 


^cc. No. A. 104133 

” • MM 


Therefore , 

In 1 yr . , 


0.326x10^ gm of Mg is in 
3 

0 0795x1 0 cm'^ of water 
total amount of water by rain fall 

6 3 

= 3.5x1x10 cm 


Therefore, time required = 0.0795x10 /3.5x10 years 

6 

= 0 . 0227x1 0 years 
= 0.0227 m y . 

The time of weathering of a 26 m thick gneissic parent rock 
and Its replacement by an equally thick layer of bauxite 
and/or latesite, calculated on the basis of chemical mass 
balance of four soluble catious (vu, Na, K, Ca and Mg) are 
presented in Table 3.8. 

The above data indicate a rather narrow time range for 
weathering Thus, in the present case, an acceptable time 
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for weathering of a 26 m thick parent rock can be around 
0,1 m.y. However, in our c a 1 cu 1 at i ons , we have considered 
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TABLE 3.8 


For element 

Ca 1 cu 1 ated T ime 


(million years) 

Na 

0 1 

K 

0 096 

Ca 

0 0256 

Mg 

0.0227 
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that all of the ram water have infiltrated in the rock. 
But, surface run off would constitute a considerable amount 
of rain water Hence, decrease in amount of water, which 
acted as leaching agent, would contribute in raising the 
time of weathering in all the above calculations. Thus, 
the times we have calculated are underestimates. 

5 EVIDENCE FROM MICROSCOPIC VIEW, KARNATAKA BAUXITE 
DEPOSIT 

In this section, we extend the calculations to 
microscopic scale. Fig. 3,3, reconstructed from the figure 
provided by Khanadali and Devaraju (1987) shows a central 
unaltered feldspar grain and the gibbsite band surrounding 
it along the periphery Here, we can assume direct 
gibbsitic alteration of plagioelase feldspar. We attempt 
to calculate the time required for the formation of this 
gibbsitic band around the feldspar grain, by means of 
chemical mass balance. 

To ensure accuracy, we have measured the thickness of 
the altered rim in six places and the average was taken. 
Average thickness of the gibbsite layer = 0.65 cm 
(measured on Fig. 3.3) 

Magnification = 70 times 

Therefore, actual thickness = 0.65/70 = 0.00928 cm 

When all other parameters remain constant, length of time 

of weathering becomes proportional to the thickness of 



UF - Unaltered feldspar 


Schematic diagram drawn from photomicrograph 
showing transformation of plagioelase feldspar 
(centre) into gibbsite along grain boundaries, 
fractures and cleavage planes {)C70) 

(After Khanadall and Oevraju, 1988) 
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TABLE 3.9 


For el emen t 

Time of weathering (months) 

Na 

4.34 

K 

4.17 

Ca 

1.11 

Mg 

1 .00 
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weathered 1 ayer . 


We 

can 

u 1 1 1 1 z e 

the 

calculations of 

the previous 

section 

for 

release 

of 

various ions from 

a 

26 m th 1 ck 

parent 

rock. 

F 0 1 lowing 

this principle. 

the 

times of 


weathering of a feldspar grain, through a thickness of 
0,00928 cms, for different ions are listed in Table 3.9, 
Geologically speaking, the above data indicate quite fast 
rate of weathering, inspite of the fact that these time 
values are underestimates The most plausible reason of 
this fast rate of weathering can be attributed to extremely 
high rainfall (3500 mm, annually) in the area from which 
the sample was taken (South Kanara district, Karnataka, 
India), 

The following schematic diagram (Fig. 3.4) shows the 
idealised picture of weathering of the feldspar grain 

The above discussions and calculations show the time 
required to weather a rock of given thickness to its 
weathered product of equal thickness under various natural 
physico-chemical conditions The variables include rock 
type, thickness of weathering, rainfall, prevailing P-T 
conditions, pH etc. 

To sum up the above findings, all the relevant data 
and the derived time values are furnished in a tabular form 
below This serves the purpose of having a comp ar at i ve 
view as well as rate of weathering in each case. 

From Table 3.10, it can be seen that the rate of 
weathering varies drastically from region to region mainly 



Altered gibbsitic layer 


UF- Unaltered plagioclase feldspar 


Fig. 3.4 



Idealised diagram of weathering of a feldspar 
gram showing time of iseathering 
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depending upon rainfall and rock types. A general survey 
reveals that the rate is more in areas of heavy rainfall 
than in areas of relatively less amount of rainfall. 

Recently, Chen et al (1988) have calculated the rate 
of weathering of an de s i t i c /andes i t i c tuffaceous rocks in 
Tatun valcanic area, N Taiwan According to them, the 
weathering rate, in the area, is indicated by the time 
required for accumulation of radioactive Be in the soil 

column They calculated the weathering rate to be around 

-1 6 
12 cm Ka or 120 m/10 years. 

Comparison of this rate of weathering, measured by a 
completely different technique, with our calculations 
suggests similarity in the order of magnitude of the rates. 
This lends support in favour of the validity of the method 
of chemical mass balance as a simple but useful tool for 
quick estimate of time of weathering. 


u 
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CHAPTER 4 


DEVELOPMENT OF A MATHEMATICAL MODEL OF WEATHERING 


4.1 INTRODUCTION 

In chapter 2, various physico-chemical factors 
governing the phenomenon of direct g i bbs i 1 1 zat i on have been 
identified. Among them, one of the most important factors, 
VIZ., pH range has been derived with precision, by the help 
of equilibrium diagrams. Now, having proved the phenomenon 
of direct g i bbs i t i z at i on of feldspars to be operative in 
nature, we, here, attempt to investigate some aspects of 
the reaction mechanism. Before stating the problem, it is 
essential to mention the background of the investigation. 
Reaction mechanisms constitute one of the most significant 
but least studied aspects of weathering processes. We 
highlight the aspect of diffusion in the weathering of 
potash feldspar to gibbsite directly, in the following 
discussion. 

4.2 STATEMENT OF THE PROBLEM 


Two principal ideas have gained attention as probable 
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procesbes of mineral zonality evolution, i e., generation 
of layers of weathered products, in a weathering profile . 

1. Weathering rate is controlled by diffusion of dissolved 
components through a layer of reaction products. In this 
process, the external reacting fluid has to traverse 
through an already formed product layer to come into 
contact with the fresh reacting solid. The reaction front 
thus proceeds inward, i.e., towards the direction of 
movement of the reacting fluid within the weathered rock. 
According to different authors, reactions of this type may 
be congruent or incongruent (Wollast, 1967, Dobrovolsky, 

1 98 7 ) 

2. The other idea emphasizes that the limiting step in 

dissolution IS not diffusion but chemical reaction, at the 
mineral-water interface, particularly at structurally contro- 
lled sites (Berner, 1978 j Holdren and Berner, 1979; Colman and. 
Dethier, 1986 , ) , 

3. A third mechanism is still on the anvil, which favours 
total dissolution of the ions from the parent mineral into 
the associated solution and thereafter insitu precipitation 
of the minerals after the ions have been recombined. This 
mechanism heavily depends on the relative mobility of 
various ionic species derived from the parent minerals 

In the present thesis, we are concerned with the 
problem of direct g i bbs 1 1 1 z a 1 1 on of feldspars and its 
possible reaction mechanism. Here, we will confine 
ourselves to the investigation of the relative importance 
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of diffusion vis-a-vis chemical reaction in the process of 
direct g 1 bb s 1 t 1 z a t 1 0 n of potash feldspar. In other words, 
we consider the alteration reaction to be diffusion 
controlled and with this assumption, try to find out the 
1 mp 1 1 c at 1 1 on s of the diffusion process. 

4.3 FORMULATION OF THE PROBLEM 

1 ) Background 

According to Wollast ( 1 967 ), direct g i bbs 1 1 1 z at i on can 
take place according to the following steps 

Step 1 ; Rapid exchange between H of the reacting 

solution (water) and K of the feldspar leading to breakage 
of feldspar structure and formation of a thin layer 
enriched in aluminium. According to Wollast (1967), this 
A1(0H)3 (amorphous) layer forms all over the feldspar 
grain. However, according to others, like Berner and 
Holdren (1979), if it at all forms, it might be confined to 
the weak zones of the crystals, such as, along twin axes or 
cleavage planes of the crystal and presumably of sub- 
electron microscopic dimension, which is difficult to 

1 

detect. No matter which mechanism is true, it does not 
rule out the importance of the presently adopted technique 
(i.e., diffusion controlled mechanism) to judge the 
relative significance of diffusion over chemical reaction. 

This step IS very fast and according to many workers 
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(e.g , Wollast 1967, Luce, 1972 and Lagace, 1965 in 
Dobrovolsky (1 98 7)) goes into completion within few minutes 
in laboratory scale experiments. 

Step 2 Once a thin layer (the so called "Protective 
layer" - of Al(OH)^ amorphous) has been formed, now, silica 
starts coming out of the feldspar structure due to 
continuous reaction with the solution (water). Aluminium, 
because of its extremely low mobility is retained. In this 
stage, the rate of release of silica is controlled by two 
factors, VIZ., rate of diffusion of these species through 
the residual layer and the rate of the chemical reaction 
between feldspar and water, at the m i n er a 1 - wa t er interface. 

Following Wollast (1967), a third mechanism can be 
included, which determines the k ao 1 i n i z a t i on from the 
gibbsite, so formed . 

Step 3 • Once the external silica concentration reaches a 
maximum value (dependent on pH), it starts reacting with 
the already formed A 1 ( 0 H )3 (amorphous) and gives rise to 
kaolinite ( A 1 2 ^ 12^5 ( 0 ) 4 ) Here, silica reacts in the form 
of silicic acid (H^i 0 4 ) 

However, for the present, we shall confine ourselves 
only within the step 2 . 

11 ) Qualitative formulation* 

In the adjoining diagram (Fig. 4 1), ABCD is the 
cross-section of a potash feldspar block After certain 
time of weathering, A'B'C‘D‘ becomes the cross-sectional 
area of the unweathered feldspar, while the annular region 
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of width X IS the weathered product. Here, we assume that 
the water-rock interaction is isotropic and it always 
proceeds from the periphery to the centre of the mineral 
grain This allows us to visualize the weathering process 
as a one dimensional transient problem, where, a reaction 
front ABCD (Fig 4 1) existing at any time t takes the new 
position A’B'C'D' at some later time t' (t' > t). The 
weathered annular region, in this case, attains a 
composition of AHOH)^ (amorphous) (amorphous gibbsite). 

In the present case, however, we extend the weathering 
reaction experienced by a single block of potash feldspar 
to a more general one, which involves a natural assemblage 
of several such feldspar grains 

The problem can be formulated with reference to figure 
4 2. PQRS represents the cross-section of a column of rock, 
initially rich in potash feldspar and having certain amount 
of porosity. Rain water flows continuously through this 
porous column with a certain range of velocity. Continuous 
weathering reaction is going on between the mildly acidic 
rainwater and potash feldspar. At any time t (t > 0), PQFE 
represents cross-section of the product layer adhering to 
the unreacted feldspar rich section EFRS. We assume that 
the reation is at equilibrium only along the line EF. 

pH value at the boundary PQ is different from that at 
the reaction boundary EF. In the latter boundary, it is 
governed by the equilibrium reaction Hence, through this 
distance we monitor concentration of H'*’ ions due to 


WATER FLOW 


I I i I i 


Schematic diagram showing cross-section of a 
partially weathered feldspar grain 


WATER FLOW 


i V ir 1 1 


Schematic diagram showing vertical section of 
a partly weathered feldspathic rock body 
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connvective diffusion. 

On the other hand, at the reaction boundary EF, 

silicic acid IS produced as a result of the reaction 
+ + 

KAIS13OQ+H +7H2O = A1(0H)3+K +3H4S1O4 


This silicic acid 

is carried out 

0 f 

the system 

by 

corrective diffusion. 

Hence, through 

the 

d 1 stance x 2 

we 


will monitor concentration of H4Si04Over known distances, 
111) Mathematical Formulation : 

The equation of convective diffusion which arises from 
mass balance across a plane perpendicular to which 
diffusion takes place, is • 



-F V. 




where, C = concentration of the diffusing species A 
A 

(In the present problem, these are H and H^SiO^) 

X = d 1 stance 
t = time 

V = velocity of percolation of ground water 
D = diffusion coefficient or diffusivity 
of the diffusing species A 


( 1 ) 


4.4 SOLUTION PROCEDURE 



One 

parti a 1 
Method " . 


of the most suitable methods 
differential equation is by 
In this case , an imp licit 


to so 1 ve the above 
“Finite Difference 
finite d i f f erence 
exp 1 1 c 1 1 form of 


scheme has been preferred over the 
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rt I I 

Vt ' 


,“!■ 'I t ' 


m 


I I'l 


f 


, iljl^ ll' 


solution in order to overcome the stability problem. 

We define the following dimensionless variables x*, t* 


and C^* as follows 


X* = 


t* = 


Ref 


'Ref 


r * = — 
C 


'Ref 


The abbreviation "Ref." stands for some reference value of 
the parameter 


d , ^ 

Now, ^ (C.* C. ) - C. ■ *- 

^ ^ ^Ref ^Ref ^ -“^Ref dt „ 

dc * 


dt 

3 t“ 


'Ref 


■^Ref dt 


dc 


\ef < 


' * ^ ^Ref dt 

Again, 


or 


Again , 


d / * 

^Ref 

^Ref 

'^A 

37~ 



dx"" 

_ ^Ref 

dc 

^Ref dx* 

' ' ■Sir 

^Ref 

dx 

d^C. 

A dc. 

A 

c. 

dc * 


dx 


= -51^ “ ■57 ( X 


Ref 


dx 


'Ref 


5 c 


'Ref 


dx 5 x 


A >> dx 


or 


2 

d C 


A _ ^Ref 


d^c 


dx^ 


'Ref 


dx* 


( 2 ) 


(3) 


(4) 
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Substituting all these expressions in equation (l), 


'Ref 


'Ref 


— T- + V. — 


Ref 


dt 


Ref 


6C 




^ = D 

* “LI 


= D 


A 


'Ref 


■Ref 


5,2 "K- 

^ '=A 

ax 


or 


dC, 


^Ref 6t^ 


V 


^Ref 




dC 


5,2 * 

^ '^A 


Ref 6 X 


2 

Multiplying both sides by 


(5) 


Ref 


'^Ref ‘^A 


dc 


dt 


V X. 


Ref 


D 


dC, 


dx 


^2 » 


" [ ' 

. I 1 ' 


I. ';i 

'I I ii 


dx 
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^2 * s * n2 * 

or, Pel — + Pe2 — ^ ^ 


dt 


6x 


dx 


where , 


Pel = 


^Ref 

■Ref'^A 


V X 


and Pe2 


IT 


Ref 


Pel and Pe2 are two dimensionless quantities known as 
Pec 1 et numbers . 

Now, while solving the above equation, we adjust the 


values of t 


‘Ref 


Ref 


D, 


, so that Pel becomes unity. 


As 


we are left with one Peclet number viz; Pe2, henceforth we 
denote Pe2 by Pe only. 


So, the new dimensionless form of the equation becomes 


be 


dc 


+ Pe 


s2 * 


at* ■•''ax* 


(6) 


Now, in an implicit finite difference scheme (Carnahan, 
Luther and Wilkes, 1969), equation (6) can be written as : 


AX 


or, 




At 

(a x)^ 


^ ‘^1-1 , n+l"^‘'i ,n+l'^‘"i+l , n+1 ^ 
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multiplying both sides by > we get, 

2 

“ ^i-l,n+l“^‘"i,n+l''^i+l,n+I 


or, a + Y C, . , , - d c 


where, 


1 , n +1 


, r _ ( 

^'^i+l,n+l(-l) ~ 
1+1, n+1 " ^ ^i,n 


1 ,n 


(7) 


a = - (Pe AX +1) 

P = A X + 2 

Y = - 1 

and d = ' ^ ' A coefficients of 

the generalized algebraic equation (7). 

In the above formulation 

1 - variable representing distance nodes 

n * variable representing time nodes 
AX - difference between two successive node points on 
the distance axis = distance interval 
At - difference between two successive node points 
on the time axis 
= time -step interval 


4,5 BOUNDARY CONDITIONS 


The boundary conditions can be ascribed with reference 
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to the following diagram (Fig. 4.3). In case of 
investigation of hydrogen ion concentration, this grid 
represents the zone of product layer of composition 
AllOH)^ (amorphous), while, in case of investigation of 
silicic acid, this grid represents a zone inside the 
unweathered rock body, below equilibrium reaction zone of 
weathering This is further illustrated in the figure 4.4. 

Boundary Conditions for Diffusion of Hydrogen ions . 

-5.7 

At boundary AB (Fig. 4 3), C(x=0, t=t) = 10 m/1 

(from average rainwater pH value) 

- 10.43 

At boundary CD (Fig 4.3), c{x=l, t=t) = 10 m/1. 

The second boundary condition arises from the equilibrium 
pH value derived from the point Q of the metastable 
boundary between K-feldspar and amorphous gibbsite (Chapter 
2, Fig. 2.6) 

Initial Condition for Diffusion of Hydrogen Ions : 

A linear initial concentration profile between the two 
boundaries has been used throughout the investigation in a 
dimensionless form. As expected, this does not affect the 
final steady state solution as has been checked with 
various arbitrary concentration profiles as initial 
cond 1 t 1 ons . 

Boundary Condition for Diffusion of Silicic Acid . 

— 5 5 

At the boundary AB, C{x=0, t=t) = 10 m/1. This is 


^'c 

Ir 


I' I 

r 



the equilibrium concentration of silicic acid furnished by 
the point Q(Chapter 2, Fig. 2 6) 

1/2 

At the boundary CD, C(x=l, t=t) = 2.K t 


TIME 


B 


D 


4 

n+1 




C 



DISTANCE 


• 3 Schematic r epr« '^ent at 1 on of finite difference 
grid 



Equilibrium 

reaction 

zone 


Downward 
percolation of 
ground water 


Fig. 4.4 


Schematic diagram showing vertical section of 
a partly a 1 te re d /un a 1 1 e r ed rock capped by a 
product 1 ayer 
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(ftssuming, parabolic rate law) 
where, t = time 

K = rate constant for release of silicic acid 
At the boundary CD, we have applied the parabolic rate law 
of Wollast (1967) 

Initial Condition for Diffusion of Silicic Acid 

Same arguments are also applicable here, as in case of 
hydrogen ion diffusion. 

In all the above cases, we have normalised all the 
variables to a range from 0 to 1 by dividing them with 
suitable reference values. 

DATA SET ♦ 


A. For diffusion of hydrogen ions . 

-7 

Percolation velocity of ground water (v) = 10 cm/sec. 
Diffusion coefficient for in aqueous medium (D^) 

= lO^cm^/sec (Garrels and Mackenzie, 1971) 

Boundary conditions for concentration • 

(Before normalisation)* 

* • the term "Normalisation" here means division of the 

parameters by a reference value in order to make it 
d imens i on 1 ess . 

C(x-0, t = t) = lo”^*^ m/1 

C(x=l, t = t) = 10 

C = 10“^*^ m/1 


I Hr 


' I 



Ref 

Boundary conditions for concentration (after normalisation) 

C(x=0, t=t) = 1 0 

“4 . 73 

C(x=l, t=t) = 10 ' (Both dimensionless) 
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B. For diffusion of silicic acid 

-7 

Percolation velocity of ground water(v) = 10 cm/sec 

^Diffusion coefficient for silicic acid (H 4 S 1 O 4 ) 

in aqueous solution (Dg) = 10^ cm^sec (Applin, 1987). 

^Different authors have reported different values for Dc , 
-5 -6 

ranging from 10 to 10 cm^sec. However, as majority 

-5 2 

of the authors have chosen its value around 10 cm /sec at 

25^0, we have chosen so. (Applin, 1987), 

Boundary conditions (before normalisation) 

C(x=0, t=t) = 10 “^’^ m/1 

1 /2 

C(x=l , t=t) = 2 K t ^ m/1 

for silica release at the 

I 

(Wollast, 1967) 


N.B, The value of K has been experimentally determined by 
Wollast (1967) for the pH values 4,6,8 and 10, In the 
present case, pH ?s 10 43, at equilibrium Thus, the K value 
corr espon d 1 n g to pH 10 has been chosen from Wollast's 
experiment, being the nearest). 

Therefore, boundary conditions (after normalisation) i 

C(x=0, t=t) - 1,0 


where, K is the rate constant 
outer boundary CO (Fig. 4 3). 

mg 


K = 10 


-4.5 


1 . (hour) 


172 


= 10 


-11.26 


moles 


1 . ( sec) 


TPl 
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1/2 - 5.5 

C(x=l , t=t) = 2 K t /10 

(Both dimensionless) 


4.6 RESULTS AND DISCUSSION 


The results were obtained by solving the appropriate 
finite difference equations in an interactive DEC-1090 
computer housed at the computer centre of IIT Kanpur. All 
the relevant diagrams were also plotted in the same system, 
using a PLOT-10 package. Programmes have been written in 
FORTRAN language and have been incorporated in the 
append i x . 3 . 

(a) Diffusion of hydrogen ions through three separate 

distances, with percolation velocity remaining 
constant 


Through three arbitrarily chosen distances (i.e., 

thicknesses of the product layer) viz- 200 cms, 400 cms and 

+ 

600 cms, diffusion of H was investigated under the given 
boundary conditions. In all the three ^cases, velocity of 
ground water was kept constant at 10 cm/sec. Table 4.1 
displays the result of diffusion of . Peclet numbers in 
the three cases and calculated time to attain steady state 
are also shown. Figure 4.5 di agr ammat ically shows the 
steady-state concentration profiles for the three cases. 

It can be easily seen that as the value of Peclet 



r;^ 

i,r 




DIMENSIONLESS CONCENTRATION 



Fig. 4,5 Plot of steady state concentration profile of 

H for three different Peclet numbers 


4 
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TABLE 4.1 


Case 

Distance 

Distance 

Pec 1 et Number 

Calculated time to 

No 

(1) 

(Reference ) 

(Pe=V 

attain s t eady- s t a t e 


( cms ) 

Ref ) 

( D 1 mens i on - 

(years ) 



{ cms ) 

less) 


1 . 

200 

200 

0 2 

3 . 8 


( = 1 2xlO*^sec ) 

2 400 400 0.4 15.22 

{=48x10^ sec ) 

3. 600 600 0.6 34.25 

(=108x10^ sec) 


u 
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numbers increases, the curvature of the steady-state 
profile also increases correspondingly 

(b) Diffusion of silicic acid through three separate 
distances, with percolation velocity remaining 
cons t an t 

As in the case of hydrogen, in this case too, three 
separate distances (i.e,, thickness of partially weathered/ 
unweathered rock) were chosen viz. 50 cms, 200 cms and 400 
cms for diffusion of H 4 S 1 O 4 under the relevant boundary 

conditions. Velocity of ground water was kept constant at 

-7 

10 cm/sec. 

Table 4.2 displays the result including Peclet number 
and time to attain steady state for the three cases. Fig. 
4 6 d 1 ag r amma 1 1 c a 1 1 y shows the steady state concentration 
profiles over normalised distances for the three cases. As 

before, it is easily discernible that with increment in 

« 

Peclet number the profiles progressively become more 
curved, 1 e., tend towards attaining the initial highest 
value (i.e,, 1 0 ) of one end. 

As expected from the theory, a high value of Peclet 
number would lead to the formation of a boundary layer. 
Although, the values of Peclet numbers used in these cases 
were not sufficient to show any concentration boundary 
layer in the real sense, the results, however, definitely 
show the real trend. In all the above cases, velocity of 


DIMENSIONLESS CONCENTRATION 


Pe = 4 0 


.Pe=20 


\P2=0-5 

s X 
N \ 

\ \ 


Initial profile-"'^ 


DIMENSIONLESS DISTANCE 


Plot of steady state concentration profile of 
H4Si04 for three different Peclet numbers 





74 




TABLE 

4 . 2 


Case 

No 

Distance 

(U 

( cms ) 

Distance 

( Re Terence ) 

(cms) 

Pec 1 e t Number 

(Pe=V.Xj^g^/Dg) 

(Dime ns lorless 

Ca 1 cu 1 ated time to 

att a I n steady state 

(years ) 

) 

1 

50 

50 

0.5 

3.17 

(=10^ sec) 

2 

200 

200 

2.0 

38.05 

(=12x10^ sec ) 

3 

400 

400 

4.0 

152.20 

(=48x10® sec) 


'1 f 
r 

h ■ 
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water percolation (v) was kept constant at 10 cm/sec. 
So, the values of Peclet numbers become function of rock 
size only, for a constant value of diffusion coefficient. 

In mass transport calculations, the final steady- 

state profile should be independent of initial 

concentration profile. To test the validity of this 

concept in the present cases, two other initial conditions 

(1 and 2) were used for the same set of boundary 

conditions It has been found that the steady state 

profiles do not change at all. These are shown in figures 

+ 

4.7 and 4.8 for diffusion of H and H 4 Si 04 respectively 
with one Peclet number for each case For other 
Peclet-Numbers, results are similar. However, in case of 
transient calculations, the number of time steps needed to 
achieve the final steady state depends on the initial 
condition. Correspondingly, the time calculated to achieve 
steady-state becomes dependent upon the initial condition. 
Thus, the times calculated for various Peclet numbers do 
not furnish the steady-state time unless the initial 
concentration profile is known. However, in the present 
cases, with the given boundary conditions, there is little 
scope for a very wide variation of initial conditions. 
Thus, the time calculated can be considered as real times. 

The fact that with increasing Peclet numbers, i.e., 
with increase in rock sizes in these cases, the time 
required to attain steady state also increases 
correspondingly is shown in Figures 4.9 and 4 10, The 


DIMENTIONLESS CONCENTRATION 



Fig. 4.7 steady state concentration profile of H 

(Peclet number =06) for two different 
initial prof lies 
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DIMENSIONLESS CONCENTRATION 
(H4S1O4) 



^-8 Steady state concentration profile of H4S1O4 

(Peclet number = 4 . 0 ) for two different 
initial profiles 


STEADY 



lOO’O 200-0 300-0 400-0 500-0 600-0 

ROCK SIZE (Cms) 



Fig. <1.9 


Plots of time required to attain steady state 
against different rock sizes, for diffusion of 
H + 
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Fig. 4.10 Plots of time required to attain steady state 
against different rock sizes, for diffusion of 
V H4S1O4 
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former shows the case of H diffusion and the latter shows 

the case of H . SiO. diffusion. In both the cases, the 

4 4 

increment is monotonic. 

(c) Effect of Change in ground water percolation velocity 
on diffusion, with distance remaining constant. 

Effects of change in ground water percolation velocity 

+ 

on the diffusion of H and H^SiO^ have been studied by 
keeping distance constant for each case and changing the 
velocity of ground water in such a manner so as to keep the 
Peclet numbers unchanged This has has been done for 
comparison with the results given in Tables 4.3 and 4.4. 

4 " 

For diffusion of H , distance is kept constant at 20 

cms and for diffusion of H^SiO^, it is kept constant at 50 

cms. From Table 4.3 it can be seen that for the three 

different Peclet numbers, viz 0.2, 0.4 and 0.6, 

percolations velocities are 10 ^ cm/sec, 2 x1 0 cm/sec and 
-6 

3x10 cm/sec respectively. However, this change fails to 
affect the calculated time to attain steady state, which in 
same (13.88 days) for all the three cases 

For the diffusion of SiO^ , table 4.4 furnishes 
similar result. In this case, for Peclet numbers 0.5, 2.0 

and 4.0, velocities of percolations were 10 cm/sec, 4x10 

“7 

cm/sec and 8x10 cm/sec. However, in all the cases, 
calculated time to attain steady state become same (7.93 
years) . 
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TABLE 4 3 

Case 

Percolation 

Pec 1 et Number 

Ca 1 cu 1 ated t ime 

No 

Velocity of 

ground water 

(Pe=V.Xf^3f/DH) 

to attain steady- 

state 


(V) 

(dimensionless) 

( Days ) 


(cm/sec ) 



1 . 

10 

0 . 2 

13.88 

2 

2x10“^ 

0 4 

13 88 

3. 

3x10“° 

0,6 

13.88 



^ * I ‘J 
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TABLE 4 4 


Case 

Percol at ion 

Pe c 1 e t Number 

Calculated time 

No . 

velocity of 

ground water 

(Pe=V.Xpgf/Ds) 

to attain ste ady- 

state 


(V) 

( cm/sec ) 

( D 1 mens i on 1 e ss ) 

(years ) 

1 . 

10 “’^ 

0 5 

7.93 

2 . 

4 x 10 "’^ 

2.0 

7 . 93 


-7 


3. 8x10 


4 0 


7.93 


I 


■) 

j. 

%'\ , 
I 

i^fi' I 

'Hr 

'hA 
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These calculations show that as long as Peclet numbers 
remain the same, any change in percolation velocity of 
ground water or distance fails to affect the time required 
to attain steady state. 

(d) Concentration profiles drawn over real distances • 

To see the dependence of steady state concentration 
profiles on rock sizes, for a given ground water velocity, 
concentrations (in dimensionless form) have been plotted 
against real (i.e., d i men s i on a 1 1 sed ) distances (i e., rock 
thicknesses) in figures 4.11 and 4 12. The former shows 
the case of diffusion of H while the latter that of 

H 4 S 1 O 4 . In both the cases, ground water velocity was kept 

-7 

constant at 10 cm/sec and the only variable was rock size. 
It can be seen that in both the diagrams, concentration 
profile IS different for different rock sizes and no part 
IS common among them. This proves that, at least for small 
thicknesses, rocks of different thickness weather in 
different ways. 

(e) Study of diffusion when both distance (x) and ground 
water velocity (v) vary. 

Figures 4.13 and 4,14 show the effects of small 

-7 

changes in the velocity of percolation from the value 10 

cm/sec on the steady state profiles, over real rock 

+ 

thicknesses. Fig. 4.13 shows the case of diffusion of H . 
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Fig. 4 11 Steady state concentration profiles of H ove 
different real distances, for three different 
Pec 1 et numbers 
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Fig. 4.12 


Steady state concentration profiles of H 1 0 ^ 
over different real distances, for three^ ^ 
differ^ ent Peclet numbers 
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For the same rock thicknesses of 200 cms, 400 cms and 600 
cms , the effect of two different velocities on steady 
states are shown It can be observed that for a small 
change in the magnitude of percolation velocity (v), the 
steady state profiles do not change appreciably for smaller 
Peclet numbers However, a large change in percolation 
velocity, with the thickness (i.e., distance) remaining 
same, will increase the value of Peclet number appreciably 
and as a consequence, steady state profile will also change 
markedly. 

Fig 4.14 illustrates the same effect, as described 
above, in case of diffusion of H 4 Si 04 for three values of 
rock thicknesses, vfz: 50 cms, 200 cms and 400 cms. 

(f) Profiles over large-scale real distances • 

Till now, we have been dealing with diffusion of H 
and H 4 S 1 O 4 through arbitrarily chosen rock sizes (i.e., dis- 
tances). Now, we apply the same boundary conditions to a 
problem of diffusion of H and H^SiO^ through a rock 
thickness of 70 m. This thickness has been chosen because 
in Chapter 3 of this thesis, we have calculated the time of 
weathering of a 70 m thick parent rock by chemical mass 
balance technique. So, it will help in comparing the 
results with the steady state time obtained by the method 
of diffusion. It may possibly suggest the relative role of 
diffusion in the total time required for the formation of a 
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water percolation velocity 
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Fig* 4 Variation in steady state concentration pro- 

files of H4S1O4 for different values of 
ground water percolation velocity 



given thickness of weathered product. 

Fig ^.15 shows the steady state profile of concentra- 
tion of against distance, i.e,, rock size (both in 
dimensionless form). Fig. 4 16 shows the same parameters 
for the diffusion of H^SiO^. Table 4.5 furnishes all the 
relevant data with the calculated time required for 

attaining steady-state. 

From table 4.5, it can be seen that although the 

distance through which diffusion takes place, is same for 
both the cases, there is a marked difference in 

steady-state time. This is because of the difference in 
the value of diffusivity. Distance, i.e., rock thicknesses 
remaining same, the less the value of diffusivity, the more 
the time required to attain steady state. 

A notable feature of these two diagrams is the 
existence of concentration boundary layer. For diffusion of 
H (Fig. 4.15), it can be seen that upto a distance of about 
(0.62x70) m, i.e,, 43.4 m, concentration at every grid 

point IS same as the upper boundary value (= 1.0, in 

dimensionless form) and while, for diffusion of H^SiO^ 
(Fig. 4.16), upto a distance (i.e., thickness) of about 
(0.88x70)m, i.e., 61.6 m, concentration is same as in the 

upper boundary (= 1.0, in dimensionless form). Beyond 

these thicknesses, boundary layers seem to exist where 
diffusion has a significant role in creating concentratiion 
gradient. In the former case, thickness of the boundary 

layer (d) is 26.6 m while in the latter case, it is 8.4 m. 


DIMENSIONLESS CONCENTRATION 



DIMENSIONLESS DISTANCE 


Fig 4,15 Diagram showing concentration boundary layer 
for diffusion of H**” (Peclet number =7 0) 
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Fig. 4 16 


Diagram showing concentration boundary layer 
for diffusion of H 4 JS 1 O 4 (Peclet number = 70.0) 
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TABLE 4 5 


,ase 

1 

Distance 

( 1 . e . , rock 

s 1 ze ) ( 1 ) 

(m) 

Distance 

reference 

(m) 

Pece 1 et Number 

(Pe=V.Xj^g^/D) 

(Dimensionless) 

C a 1 cu 1 ated 1 1 me 

to atta 1 n steady 

state 

(years ) 

Diffusion 

of H"*" 

70 

70 

7 0 

7769.0 

Diffusion 

of 

H 4 S 1 O 4 

70 

70 

70.0 

77689.0 



I 




“■r 

i I ‘ 

W 'f ' < 1 
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4 7 COMPARISON OF NUMERICAL SOLUTION WITH ANALYTICAL 
SOLUTION OF THE GOVERNING DIFFERENTIAL EQUATION 


The differential equation for convective diffusion is: 

V. ^ LI' — O (8) 


+ V. ^ = n D- ^ 
^ dx2 


where, c = concentration 
t = time 
X = distance 
D' = diffusivity 
t| = tortuosity factor 
V = percolation velocity of ground water 


In our present study, 
a separate parameter, 
value of an effective 
account the effects 
equation changes to • 

+ V. - D 


we have not considered tortuosity as 
On the contrary, we have used the 
diffusivity (D) which takes into 
of tortuosity. So, this governing 





(9) 


Golubev and Garibyants (1971) have discussed about the 
analytical solution of the above equation. The following 
account has been largely drawn upto from their treatment. 
They have treated the problem of convective diffusion 
through porous media with various boundary condition. The 
most appropriate of them, for the present study is the case 
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of - "Filtration from a steady source" This is because, 
in our cases, we have always maintained the concentration 
at one end fixed (i e., c(x=0, t=t) values are constants 
for the diffusion of h"*" and H^Si 04 ) 

Thus, under the boundary and initial conditions of* 

X >0, t=0, c(x,t) =0 
X = 0, t > 0, c(x,t) = constant 

= (say) 

The solution of the equation (9) has the form : 

C(x,t) = ^ L( 1-erf (X=1LL.)) + ( l_erf )] (10) 

2fDl 2f^ 


where , 


erfx IS the error function of x defined as 


erfx = 


7 ^ 


/ exp(-y^) 


dy 


In the subsequent analytical treatment, we take the initial 
concentration profile as given in the condition above, 
i.e,, for x > 0, t = 0, C(x,t) = 0. 


Comparison between analytical and numerical solutions - 
for diffusion of hydrogen ions (H*^ ) : 



Out of the three distances (i.e., rock sizes) which we 
considered in numerical treatment, only two values, viz., 
200 cms and 600 cms have been chosen for which diffusion of 
.H"*" 1 s compared . 
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For the present cases . 

+ 2 

D = diffusivity of H = 10 cm /sec. 

t = time required to attain steady state for diffusion 

through rocks of 200 cms and 600 cms length 

-7 

V = 1 0 cm/sec . 

From the numerical result (obtained with the help of finite 
difference procedure), time required to attain steady state 
IS first computed. This value of t along with the values 
of D and V are substituted in the equation (10). Defining 
dimensionless concentration variable C* as C* = C(x,t)/CQ, 
we arrive at a functional relationship between C* and the 
distance variable x. Figures 4 17 and 4.18 display the 
plots of C* against x (in cms), for two rock sizes viz* 200 
cms and 600 cms respectively. In each plot, the dashed 
line shows concentration profile drawn from analytical 
result while the solid line shows concentration profile 
drawn from numerical results. 

It can be noted that in both the cases, trends of 
concentration distribution are- similar and comparable. 
However, towards the boundary x = 1 (= 200 cms and 600 cms 
respectively), increasing disparity exists between the two 
profiles. This disparity is expected and can be explained 
in terms of dissimilar conditions at the boundary x = 1. 
In case of analytical treatment, solution was essentially 
obtained for a semi-infinite solid with no condition 
specified at x = 1, while, in case of numerical method, 

there is always a fixed concentration of H specified by 
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the prevailing equilibrium reaction at x = 1 . 

Comparison between analytical and numerical solution - 
for diffusion of silicic acid 

Out of the three rock-sizes which were taken into 
account for numerical investigation of diffusion of silicic 
acid, two cases, with rock sizes of 50 cms and 400 cms are 
now considered for comparison with analytical solution. 

The same initial condition as taken before is 

applicable here too as well. 

For these cases . 

0 = diffusivity of H^SiO^ = 10 %m2/sec. 

t = time required to attain steady-state 

-7 

V = 10 cm/sec. 

As done before, at first, time required to attain steady- 

state for diffusion of H^SiO^ is computed, numerically. 

With known values of t, V and 0, we get a functional 

relationship between C* and x (in cms) in these cases also. 

Figures 4.19 and 4.20 show plots of C* against x for 

diffusion of H^SiO^ at steady-state time. 

+ 

As in case of H , here also a progressive, disparity 
between the numerical and analytical profiles is observed 
as we approach towards the boundary x = I (=50 cms and 400 
cms respectively). In these cases, the reason of this 
disparity is the presence of a time dependent condition 
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Fig. 4.20 Comparison of numerical and analytical results 
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(Distance taken in numerical solution=400 cms 
Pecklet number = 4 0 ) 
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{ Wo 1 1 ast-r ate of release of silicic acid) at the boundary x 
= 1, which IS ignored in case of analytical solution and 
once again sem i - 1 nf i n i te formulation is used there with no 
specified condition at x = 1 

Calculation of maximuni distance traversed by the diffusing 
substances under semi-infinite mode of diffusion • 


With some reasonable assumptions, equation (10) can be 
expressed in such a way so that the maximum distance that a 
diffusing species will travel, can be computed The 
following treatment is adapted from Golubev and Garibyants 

and t are rather large, then expanding 


(1971). If X 
X + Vt 


erf ( 


2vm 


) into a senes, it can be shown that the 


second term of the equation (10) is small as compared to 
the first term and may be neglected Therefore, the 
asympto-tic solution of equation (9), for conditions as 
described before, has the following form 


C( x,t) = — S [l - erf ( 


X - Vt 

2fm, 


)] 


(11) 


Without any great error eqn (11) can be written as , 

C. 


C( x,t) = ^ [l - 


X - Vt 

YlTDt 


] 


( 12 ) 


Now, for a minimal concentration which is still detectable 
in practice (i.e., chemically), we have, 


"min 


= [1 - ( 


max 


- Vt 


■yTciSt 


-3 


(13) 
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Among the two parts of the expression for ^ max ’ 
corresponds to the d i s tance -component contributed by 
percolation velocity of the solution. The other part i.e., 
V^Dt corresponds to the distance-component contributed by 
d 1 f f u s 1 on only 

The parameter t in the expressiion of j can be 

any time after the onset of diffusion. However, in present 
cases, we have used only the time required to attain steady 
state in each case. Using these time values along with 
appropriate values of V and D, theoretical values of x 

max 

are calculated. 

Tables 4.6 and 4.7 furnish the distance of diffusion 



(i.e., rock sizes) as taken in numerical analyses along 
with the calculated values of x 

max 

Figures 4.21 to 4.24 show the concentration 
distribution upto theoretical x along with the 

concentration distribution furnished by numerical solution 
upto the pre-defined distances. * The first two plots show 
the case of diffusion of H while the last two show that 
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TABLE A. 6 


Diffusing substance Hydrogen ions ( H ‘^) 


Case 

Pec 1 et Number 

D 1 s t anc e ( 1 ) 

Calculated 

max 

No 

( D 1 mens 1 on 1 es s ) 

(cm) 

using eqn (15) 

( cms ) 

1 . 

0.2 

200 

372 

2, 

0.6 

600 

1422 


'I 


y, 
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TABLE 4.7 



Diffusing substance 

• Silicic acid 

(H 4 S 1 O 4 ) 

Case 

No. 

Peclet Number 

( D 1 men s i on 1 es s ) 

D 1 s t ance ( ] ) 

(cm) 

Calculated X 

using eqn .(15) 

( cms ) 

1 . 

0,5 

50 

113.5 

2 . 

4,0 

400 

2308.0 













of 
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From these diagrams as well as from earlier plots, the 
following observations can be made : 

(a) We have always investigated the aspect of diffusion 
through finite rock-size. The pre-assigned boundary 
conditions confine our steady state concentration profiles 
within the given distance (length of rock body). However, 
in case of analytical treatment, no condition is imposed at 
the other boundary towards the direction of flow. This 
essentially makes the problem semi-infinite type. 

(b) It IS noticeable that near the boundary having maximum 
concentration (= 1.0 in dimensionless forms), the 
concentration profiles displayed by both the methods are 
close to each other. This is because, at this boundary, 
conditions are same for both the methods. However, in the 
other boundary disparity increases rapidly This is 
because of dissimilar boundary conditions. 

(c) It is noteworthy that in all the cases, for the same 
time (time required to attain steady state in case of 
numerical solution) always larger than the 
rock-size through which numerical solution was obtained. 
This IS expected. Because, instead of imposing a boundary 
condition at other end (x = 1), if we would have allowed 


, 1 

the substances 

to 

diffuse 

freely 

through 

the 

rock 

, then 

tor 1 

1 , 1 

^ II 

y 1 1 

this X 

max 

wou 1 d 

have 

been 

such a 

distance 

where 


concentr at 1 on 

1 s 

sma 1 lest 

as 

comp ar ed 

to 

the 

source 


,i{ 



strength, but still measurable. 
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(d) The analytical solution, however, furnishes the actual 
distance through a rockbody upto which diffusion is 
significant in case of convective diffusion from a steady 
source along with a free boundary. In numerical solutions 
between fixed rock-sizes, diffusion has significance 
throughout that rock body before the formation of a 
concent r a 1 1 on boundary layer Once a boundary layer is 
formed, after that diffusion is only significant within 
that layer in generating a perceptible concentration 
prof 1 1 e . 

(e) In reality, the analytical solution with sem i - i nf i n i t e 
mode of conditions corresponds to an open system, where 
reaction products (e.g., silicic acid in the present case) 
are continually lost in the surrounding medium. On the 
other hand, numerical solution with fixed boundaries is 
more appropriate in case of a closed system (e.g,, 
laboratory simulation of weathering of a mineral grain) 
where rate of release of products can be measured at a 
fixed a boundary. In natural system, such as ours, as the 
ground water percolation velocity is very slow, for all 
practical purposes, within short period of time 
(geologically) numerical technique with fixed boundaries 
are applicable. This is because within short period of 
time not much of material will be out of the system. This 
IS further favoured by the very slow movement of the 


reaction front. 
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4.8 GEOLOGICAL IMPLICATION OF THE RESULTS 

We started our investigation in order to find out the 
role played by convective diffusion in the process of 
weathering. On the way, we have picked up a specific 
example, viz; direct gibbsitic alteration of feldspars. 
All the boundary conditions considered are consistent 
within the thermodynamic framework of direct g i bbs i t i z at i on 
of potash feldspar. Thus, the results obtained highlight 
the diffusive aspect of the transformation. Here, we have 
not considered the aspect of chemical reaction rate 
Conclusion regarding this will however appear logically. 

From the analyses that have been carried out so far, 
the following conclusions emerge out ■ 

1. Assuming equilibrium reaction between rock and water, 
precise measurements of parameters like concentration, pH 
etc can be made of the diffusing species. 

2 . For a given thickness, when steady-state will be 
attained, it appears, that the reaction will come to a 
standstill and no further weathering will take place. This 


srt , 

SJ, 1 / 

fj 
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contention, however, is not true In our treatment, v/e have 
not coupled the problem of boundary movement. But, >n 
reality, the reaction boundary will be moving at a definite 
rate (however small the rate may be). For this reason of 
continuous change in distance of effective diffusion, 
reaction will never cease, until whole of the rock is 
weathered or the weathering agent ceases to operate. 
Another reason for cessation of the weathering reaction may 
be the formation of an optimum thickness of the product 
layer which, due to its increasing resistance, hinders 
further diffusion significantly 

3. As we have outlined, at the very beginning, two aspects 
play significant role in the generation of weathered 
profiles. These are - (a) reaction rate and (b) diffusion 
of reacting species as well as of reaction products, 
through the rock mass. In the present case, for distances 
around few hundreds of centimetres, the steady-state time 
values are of the order of few years to hundreds of years. 
On the otherhand, in natural weathering of feldspars, the 
time required to weather through few hundreds of 
centimetres is much more. This leads us to conclude that 
in a natural weathering process, reaction rate is slower 
than the rate of convective diffusion. It thus acts as 
rate limiting step and actually determines the time 
required for weathering of a given thickness of rock and 
its replacement by the weathered product of equal 
th 1 ckness . 


CHAPTER 5 
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CONCLUDING REMARKS 

The present work has dealt with the process of 
weathering of feldspars with particular reference to its 
direct conversion to gibbsite. Three aspects of feldspar 
weathering have been investigated. These are : 

(a) Cry s t a 1 1 0 - chem 1 c a 1 changes and thermodynami c viability 
of weathering of feldspars to gibbsite, both directly and 
1 nd 1 rec t 1 y . 

(b) Calculation of time of weathering by chemical mass 
balance of four soluble elements (Na, K, Ca and Mg) from 
their fractionation in the parent rock, weathered product 
and ground water. 

and, (c) Role of convective diffusion and chemical reaction 
in direct g i bb s i t i z a t i o n of feldspars. 

Following are the conclusions drawn from the investigation. 

1, Weathering of feldspars, when viewed as a successive 
structural simplification should follow the normal course 
of, feldspar mica kaolinite gibbsite. Direct 
gibbsiti2ation bypasses the intermediate stages. 
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2. The thermodynamic stability of K-feldspar shows that 
direct conversion of K-feldspar to gibbsite is not 
possible. However, a metastable boundary has been shown to 
exist between K-feldspar and gibbsite. This implies that 
if gibbsite forms directly from K-feldspar, it will be in 
metastable f orm . 

From stability diagrams involving Na- and Ca-feldspars, 
direct conversion of these feldspars to gibbsite appears to 
be thermodynamically viable. 

3, In all the above mentioned cases, direct g i bbs 1 1 1 z at i o n 
should take place at high pH conditions (pH > 10). 

4 The rates of weathering leading to formation of bauxite 


vary widely. Mass balance calculations 

have 

shown 

that 

among different rocks 

for the 

same 

element and 

for 

different elements in the 

same rock. 

weathering 

rates 

are 


d 1 f f eren t . 

5. On an average, in tropical regions, rate of weathering 
seems to be a few metres to a few hundred metres per 
million years . 

6 . When the velocity of the percolating ground water re- 
mains constant, the time taken to reach steady state for 
any diffusing species {H**’or H 4 S 1 O 4 ) increases with increas- 
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ing rock-size. 

7. If the rock-length, through which diffusion occurs is 
kept constant but percolation velocity of ground water is 
varied within realistic range, it does not affect the time 
required to attain steady state. 

This is true for the diffusion of both H and H^SiO^ . 
For example, for diffusion of H through a distance of 20 

cms, under three different percolation velocities of ground 

-6 ”6 -5 

water (10 cm/sec, 2x10 cm/sec and 3x10 cm/sec), the 

time required to attain steady state has been found to be 
13.88 days irrespective of ground water percolation 
velocity. 

8. Plots of steady state concentration profiles over 
rock-lengths (through which diffusion occurred) show that 
the steady state profiles are different for different rock- 
lengths and no part of the profiles is common between them. 
This shows that rocks of different size weather in a 
d 1 f f eren t way . 

9. For very large rock-size, which contributes in 
Increasing the value of Pec 1 et - numbe r significantly, a 
prominent concentration boundary 1 ayer is formed. Thus, 
the zone of effective diffusion is reduced and is 
henceforth confined only within that boundary layer. In 
other portion, diffusion does not produce any additional 
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change in concentration profile as there is no appreciable 
concentration gradient. 

10 . It seems that diffusion process attains steady state 

well before the time required to weather a rock of given 
thickness This indicates that diffusion plays an 

important role only during the initial stages of the 
weather i ng process . 

11 . After the attainment of the steady state concentration 
profile through diffusion, the process is likely to be 
controlled by chemical reaction. 

Two chemical reactions are dominant in the weathering 
of K-feldspar to gibbsite. These are : 

j. + -|. 

(i) K ^ H (Surficial ionic exchange between H of 
water and K of K-feldspar) 

+ + 

(n) KAIS13O3 + H + 7H2O = AUOHlg + K + 3H4S1O4 

amorphous 

As experimental evidences suggest that the first reaction 
IS very fast, therefore the second reaction which produces 
H4S1O4 IS most likely to be rate controlling. 

SCOPE OF FURTHER WORK 


The theoretical approach outlined in this thesis can 
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be further extended and supplemented in several ways : 

1. In this work, we have assigned the initial 
concentration profile arbitrarily. However, from a 
partially weathered rock-profile (or in small scale from a 
partially weathered piece of mineral) the exact 
concentration of different species, at various points can 
be measured by using electron microprobe analyser. 

The calculated concentration profiles for a given size 
of rock can also be measured and verified by this 
instrument. All these refinement will enhance the accuracy 
of the mode 1 . 

2. In the present treatment, we have assumed regular 
geometry for the shape of the cross-section of the 
weathered zone. In reality, this may not be true. Thus, 
for irregular shapes, advanced method like FEW (Finite 
Element Method) can be used for better results. 

3. Differential equation containing both diffusive and 
chemical reaction rate terms should be developed and 
solved. This will help in establishing mathematical 
relationships between diffusion parameters and chemical 
reaction rate terms in’ case of weathering. 

4. Experimental work should be carried out to calculate 
rate constants for important weathering reactions. 
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APPENDIX 1 

List of minerals, their composition and Ag 2 values : 

Mineral Composition G° in k.cal/mole 

( 1978)* 


1 • Quartz 

S 1 O 2 

-204.853 

2. Gibbsite (crystalline 

)A1(0H)3 

-276.289 

3 • Gibbsite( Amorphous ) 

A1(0H)3 

-271.9^ 

4. Kaolin ite( crys tal line )Al 2 Si 205 ( 0^)4 

-908.939 

5. Kaolinite (Poorly 

crystalline) 

Al2Si2O5(0H)4 

-901.0^ 

6 * K-mica (muscovite) 

KAl 3 Si 30 j^Q( 0 H )2 

-1339.873 

7, K-feldspar (micro- 

clinne) 

KAlSi 30 g 

-895.294 

8 . Albite 

Isja AlS i^Qg 

-887.972 

9. Anorthite 

CaAl2S ^ 2^8 

-961.068 

10 . H 4 Si 04 


-312.918 

11 . H 2 O 


- 56.732 

12 . 


- 67.581 

13. Na'*' 


- 62.656 

14. Ca'*’'*’ 


-132.426 

15. H"'" 


0 

16. A1(0H)4 


-310.2 

17. Al^'*’ 


-116.0 

* - Data from Robie et 
into k.cal/mole by 

al.(1978), recalculated 
p.R. Bhavana (1983). 

from j/moie 


1 - From Krauskopf (1979) and the references therein 
^ - From Kittrick (1966). 
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Abbreviation ; X/n = Crystalline ; Amorp:= Amorphous 

1. Kaolinite (X/n) - Gibbsite (x/n) equilibrium 
Al2Si20^(0H)4 + 5H2O = 2A1(0H)3 + SH^SiO^ 

log = -10.4 

2. K-mica - Kaolinite (X/n) equilibrium 

2 KAl 3 Si 30 j^Q( 0 H )2 + 2H‘*' * 3 H 2 O = 3Al2Si205(0H)4 + 2K'*' 

log K2 = 8.82 

3. K-mica - Poorly crystalline Kaolinite equilibrium 

2 KA 13 Si 30 j^q( 0 H )2 + SH'*’ + 3 H 2 O = 3Al2Si205(0H )4 + 2K‘‘’ 

log K3 = -8.636 

4. K-mica - Gibbsite (X/n) Equilibrium 
KAl3Si30j^Q(0H)2 + H'*’ + 9H2O = 3 A 1 ( 0 H )3 + 3H4S1O4 + 

log K^ = -11.187 

5. K-f eldspar-K-mica Equilibrium 

3 KAlSi 30 g + 2H‘^ + I2H2O = KA13 Si30j^q(0H)2 + 6H4S1O4 + 2K.‘^ 

log = -10.354 

6. K-feldspar - Kaolinite (x/n) : Equilibrium) 

2KAlSi3O0 + 2H + 9H2O ” Al28i205C^H)4 ^K + 4H^SiO^ 

log K^ = -3.961 

7. Kxfeldspar-Poorly Crystalline Kaolinite Equilibrium 

2KAISI3OQ + + 9H2O = Al2Si205(0H)4 + 2K'‘' + 4H451O4 

8. K-feldspar - Gibbsite (Amorph.) Equilibrium 

KAlSigOg + + 7H2O = A 1 ( 0 H )3 + + 3H4SIO4 

log Kg = -10.398 
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9 «> K“f elds par - Gibb site ( X/n ) Equilibrium 

KAlSi30g + H'*' + 7H2 O = A1(0H)3 + k'*' + SH^SiO^ 

log K <5 = 7,18 

10. Poorly Crystalline Kaolinite - Gibbsite (Amorph.) Equilibrium 
Al2Sl2O^(0H)^ + 5 H 2 O = 2A1(0H)3 + 2H4SiO^ 

log -11.014 

11. Albite-Kaolinite (X/n) Equilibrium 

2NaAlSi30g + 2H'*' + 9H2O = Al2Si205(0H)4 + 2Na'*' + 4H4S1O4 

log X = -0.44 

12. Albite-Gibbsite (X/n) Equilibrium 

NaAlSi30g + H'*’ + 7H2O = A1(0H)3 + Na'*' + SH^SiO^ 

log K2^2 = "5.423 

13. Anorthite-Kaolinite (X/n) Equilibrium 
CaAl2Si20g + 2 h'‘‘ + H2O = Al2Si205( OH >4 + Ca'^"'' 

log K j_3 = 17.27 

14. Anorthite-Gibbsite (X/n) Equilibrium 

CaAl2Si20g + + 6H2O = 2A1(0H)3 + Ca’*"'’ + 2H4Si04 

log = 6.87 

N.B; K-feldspar and K~mica (muscovite) are always considered 
to be in crystalline state and so also albite and 
anorthite. 
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APPENDIX 2 


Dissolution of various ionic species of aluminium 

Reaction in acid range ~ 

1. A1(0H)3 + SH'*’ = + 3 H 2 O 

(crystalline logKj^ = 7.263 

g 1 bbsite) 

Reaction in alkaline range - 

2 . A 1 ( 0 H )3 + H 2 O = A 1 ( 0 H )4 + H'*' 

(crystalline log K 2 = -16.731 
g ibbsite) 


* Data taken from Appendix 1 
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nME=ri4c+3'r'r42' 

TLIM = ILtH.U 
00 30 I =3»i*G 
A(I)sAL3MV/ttorV 
BjiiiBErAyO'iAviri'’ 

:c{I)=gamaa/)cl'ia 

Af2)s0.3 
B(?) saer'AvocoiTA' 

CC(2) = GA4'v|A/35l.'r* 

A(H-l)=ALp'IA/D:LirA 
B(N-I ) = 3:TA/3CSirA 

'•»' ( ij - 1 ) s j ,3 

5TM}s2.o>5<<‘SORr(riME)/:Br 
C(I )=1,3 
C(M)=0,0 
C(M)=0. 106E-I1 

00 50 I s 3,^-2 
D(T)cC(I) 

0(2) =5 C 20 -1(A)jP HA'/ delta) •'C ( 1 ) 

D(M-i)=:i('j“i )-(gahma/oflta)*:(N) 

CALL TRIOaGC 2;, M-a , A,B,c:i, 0,0 
TrctiLiM/'isrc)^) ♦isr&p,N2 - iLim) go to 40 
WRirE(22,2D3) (:i(I),I=1,N) 

IF (TIMO’.tjE.'l ,5) G.3 ID 10 
rfRrref24;2333 

rfRire(24 

FDRMAtIBHO r'MCI, 13R,4 

1 TS) 

FORMATJH ,:iU.6//ClH , 7X , 21E1 7 .6 j ) 

foRMAr(n>cf^ aIli'^hJ^»13X, 'BETA', 12 X, 'GAHMA'fl 3x, 'A( D 'B(I) ' , 

STOP 

plMEMSI3M^A(20)}1(2jjl’i:i 25);6lio! ,V(20) fBETdOl ) ,GAM{ 131) 

)/3G.r( JF) 


,A(3),B(3),CCC3) 

ES OF n&NCEMTRAflDN AT ALL GRlO-pjIM 


BET(IF)=3(I«') 
GAM(IF)=3'( IF 
IFPlPlFtt 
00 1 I =lT?lfL' ^ 
BET(I)=dh):Al(r)t"C 
gam(i)=(: 3 i(i)-ft(n*G 


C(I-1 )/B£T(l-l) 
3i<t) -AintGAMd-n ))/9E1(I) 
V(L3=GA4 (lO 
LASlsL-Is- 
DO 2 K =1,L'A15T 

V{Sf) =3AMI( I ) -CCK I ) t=V t in ) /BEK I) 
RETURN 
END 



TriT*; p'»::;i‘-'5' PLOTS Di''E'J5nN[,L55 cnNfENrRfcmw ^gim'jst 
oi'Ie^isidnlcss TrsrA\c: 

PRDat«A''’i:T ^^RT l\ PRAri'5 OS 

DIvis^I5I3v XC25),yC 25,25i 

REAO(2i,») 

REAp(24, M ( (i (J,, = I ,21 ) , 

N = 21 ;XSlElP = 5,Cl5;y5TEP=/) 05 
X( 1 )=0,3 
00 20 1? = 1,20 
X[TP+l)sX(I?)f0.05 
fO.'l 

X^T'rsxNhO.'i , 

:ALI, GRSTRTCiDlO,!? 

Uoz 121? 

FOR'-iAi CH, 'W'liUL' c:hs:k » 0 ; acclpt +,i 

2Ai,b rii oo'»C-K'Ji>i^5.o,x^in,-rvJjrj/5,o,yJiiv) 

CALL VW3JRrC-K^l'I/'5,0,XA'IV,-irrtlN/5.0,Y/JIU) 

CALL MO/E'(0.'),0.0) 

Call oraw(Min,3.0) 

CALL, DRAWCX'J IN, IfJll'O 
CALLlORA^n.OjUl'O 
CALL DRA/iiO.O , O.O) 

CALL SMOOTH 
00 5 J = !,<’ 

CALL SKIP 

CALL DRAW(1(?^0 ,y(i , *1) ) 
coMriNu: 

CD*jriKu: 

CALL 'ID/CCO .0 , 0,0) 

CALL ve:r-:li 

yi=yrfiN?ioo.o 

DO 25 1 = 1 , 5 

CALL MOO(4.'5*X5TSP,0.0) 

CALL DRA^(0.O,-iyn 
CALL 

CALL TXrCU^n hCALL TXPCURO) 

AsXSTEP^ELpOAind^A) 

CALL R‘iJ'l^‘lC-it,l,4) 

CALL M3i/:(0.O,5.Ofyi) 
coMriTjii: 

CALL 'lO/Cf-l ,O,0.-3J 

Xl=X^/IM/lDO.'i 

DO 30 J s 1,5 

CALL HO/:(0,'3,4.OtirSi:EP) 

CALL DRA<(-Xlj,0.0) 

CALL »40V: ■^.fi^itlfO.'O) , 

CALL TXlC'JRO) f :ALti TXFCURO) 

BsysrEPdLOAif'I JiM) 

CALL R’4J-<8RU',1 ,4) ^ 

CALL H0;2(3.'>» Xl , 3. 0) 
coNriNu: 

LiIj VF*'^^S 

?AlL '4Di?:?X*l[W/2.0, -5f/IT V5.0) 

CALL TEXTl/Mt^O'I'ICNSnilLeSS 0I51ARCE') 


CALL rEXr'(l3, 'O'INJCNSIONLESS TIME") 

CALL MO/Ci , YftI')/2.03 
"ALL TXAi'JGLpCSO.'I) 

CALL TEXrn/f ’Ol'lC^SlORLeSS COMCENTRATTON') 
CALL CiRSTOP; ACCF^J ♦, J‘ 

STOP 

END 



1 JO 

no 

I'ZI 


2JO 


300 

c 


400 


5J3 


$00 


THIS PRD33V-1 PLOTS 0 I 'lF V'l DULESS CONCFV fPAT HfJ 

ASAi'Jsr on:vsronLi5En oistajce 
PRD 3q^'H:^^?\Hr^A pkau^ o: 

ST.rn rH:st5 :'OR</r,i,r aanour/i & fl7-ap 

2I'i:\)5I3^ XC25),XI(25)n(i^F25) 

Rc.flDC 22 ,n Uy C J, *13 , vi = I 21) J=1 ,3 0 
REAOC/lin).((Uf'»).*i = l/n) 

N = B{ XST-,5^10 J.-'; ySTEP^O 05; K = 3, = : 

XC13"9»0 
DO 100 I s 1,7 
xciH3sy(i3KsT::p 
XAiNsxmt'i.O' 


= 21 


u, ’R:'AOy')MCCF3,r ♦,! 

nci'it-x^?i'i/5.5,x<r^,-y.^Hfj/5,o, y<’ifn 

?DRr’(-X^£')/'5.3,y^l'J,-yAlh/5 0,y,^L^) 
i/EC0»'0f 0,03 


yrfrv=yci. 13 to. 1 

XI (10=0 0 

03 no 13 = 1,20 
XI ClPtl)=Xl (IP) a, 05 
CALL, GPSiqrtmO,!) 

TYPE 121 

CALL ViVPDRn 
CALL MDi/CCO. 

CALL DRiiVIXVINfO.O) 

CALL 0RAV(X*J[N,r,;i*IJ 
CALL DRAWiJ.'O, r/JTM) 

Call DRArfco.’0,o.o) 

CALL yrc'^s 
CALL tooth 
OD 200^ JiJ = 1 , vg 
XmXKJJ.) *2)5.0 
CALL DRArfCXr'l.ytl.JjO) 

CDNTI'JIJC 
CALL SKI? 

CALL MO/:'(0. '0,3.0) 

CALL SMOOTH 

DO 300 KK = I , I'J 

XT2 = Xl(X<’)000.0 

CALL DRAXCXl'2', y(2,KK) 3 

CO^lTTVli: 

CALL SKI? 

CALL H0/:(0 '0, 0,0) 

CALL S'^OOTH 
DO 430 LL' = 1 , HH 
Xl3sXl(Li[,')*50O.O 
CALL DR^|^CX^'3', y(3,LL0 ) 

CDNTI 7'IC 

CALL Rn>^:(0.'3 , 5 .0) 

CALL YFCRCLi 

XXsX'llN/l 30,0 

DO 509 I = 1,5 

CALL »<D/S(0.'5,4.O»ySTEP) 

CALL 0RAItf(-<X,3.0) 

*ALL «DY;(-3,6tXX,t>*p) „ 

call rXIC'jRli) t CAtt, TXFCJR(3) 


3) 


B = ySTLPf-’LaA)l’i{Xf4 
CALL RXJX3R(3’,l,4j 
CALL rtO/C(9.0*XX,5. 

CoUTI^Jli:' 

CALL MOnCO.0,-11 .0) 

?Y = X^IHyi3355’'.3 
00 60D JI = 1,6 
CALL "HQY2'(i5SriEP,0.O) 

CALL pRAvh.O,-Yy) ^ 

Call mdye(5,o, - n+yy) 

CALL ixr:'j=l(3i * -WL'L TXFCJPO) 
AsXsrEPf*“LOAir( 7) 

Call rnjh9R(m, 1,5) , ^ 

CALL klO/E'(3.0,5,0*yy) 

Call Mo5o?X'(iN/?.o,-rrfi'i/$,o) 

CALL rEXTflSf^O'ISrAHCE (IN CHS)') 

CALL Mn5:(-X»JlV/6,0,yNlV2,01 

CaL'L' }Ejr?27?^5?M:,>IsrONLESS COfCFNTRATION': 

CALu GRSro®; ACiCSPT ♦, J 

910? 

END 
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